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for rotorcraft vibration reduction. The frequency-dona in formulation is 
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SECTION 1 
INTRODUCTIO N 


1.1 BACKGROUND 

One of the major difficulties in the acceptance of rotorcraft for both 
commercial and military applications is the mechanical vibration level. 

There are several significant vibration sources in rotorcraft, including 
rotors, transmission systems, and engines. Rotor and fuselage flexibilities 
amplify the effects of the vibration sources. The level and spectrum of 
vibration varies along the fuselage. The transmission of rotor blade cyclic 
ljads into the fuselage represents a major source of vibration. It is 
easily shown that the major components of transmitted forces are at N/rev, 
2N/r'ev» 3N/rev ... etc., where N is the number of rotor blades [1,21. 

Without control, the vibration levels in rotorcraft may exceed 0.25 g — 
significantly higher than those for fixed-wing aircraft. For a four-bladed 
rotor turning at 240 rpm, the vibration would occur at 16Hz, 32Hz, 48Bz, ... 
etc. At these frequencies and levels, the vibration degrades both the 
pilot/passenger ride quality and increases maintenance requirements. 
Reduction of rotor-induccd vibration is of major importance in future 
rotorcraft developments. 

For many years, mass-spring-damper systems have been used for 
rotorcraft vibration reduction. These systems are heavy and have a limited 
effectiveness range. Considerable research has recently been done on 
'active' techniques, which reduce vibration by directly changing the blade 
pitch which produces aerodynamic forces to counterbalance the vibration 
loads. This report studies two approaches that have been proposed f«r the 
active control of helicopter vibration. 
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1.2 ACTIVE CONTROL OF ROTORCRAFT VIBRATION 


Two approaches have been proposed to reduce rotor- induced vibration in 
helicopters through active control of rotor blade pitch angle — one based on 
frequency-domain analysis and the other based on a time-domain model. The 
frequency-domain approach using multicyclic feedback control is the more 
established approach. It has been studied theoretically and validated in 
certain wind tunnel experiments [3]. The frequency-domain approach uses a 
model which relates amplitudes and phases of vibration and multicyclic 
inputs. The time-domain approach uses a model which relates the time 
history response of vibration with input time history. It has been possible 
to use this model because of recent work in frequency-shaped modern control 
design methodology [4]. It has been tested on a detailed rotor systems 
research aircraft (RSRA) simulation, which includes fuselage flexibilities 
but does not include rotor aeroelastic effects [5]. Both approaches have 
been studied in this report to evaluate potential advantages and 
disadvantages. 


1.3 SUMMARY OF KEY WORK 


The key results from this research are as follows: 

(i) There arc many similarities but important differences between 
the frequency-domain and time-domain formulations of the active 
rotorcraft vibration suppression problem (see Sections 2 and 3). 
Both approaches are similar to the extent that they are 
formulated to control vibration at discrete frequencies. The 
major differences in their behaviors arise because the time- 
domain approach updates the control at each sample point while 
the frequency domain approach updates the control after several 
vibration cycles. 

(ii) System identification techniques are available for both the 
time-domain and the frequency-domain approaches. These 
techniques can be used off-line or in real-time. The system 
ideutif icati on approach used to identify the local frequency- 
domain model can be extended for improved accuracy. 

(iii) Both the time-domain and the frequency-domain approaches 
have been analyzed for transient behavior, robustness, 
susceptibility to various noise sources, and implementation 
complexity. 

2 








'C 


,,r ’ •* 

..be rm .„ *“» ~ 

<v ’ Sissr/.'s.fsssi" r*" *••■ >““*■■> -»■«- ..... 

(Vi) Further work is necessary to develop a ful ly-adapti ve time- 


1.4 SLMMARY OF RETORT 

The report is organized as follows. 

~ C t i0n - 2 ° f th ° rep0rt snmma rizes the frequency-domain and the time- 
donaxn methods and their variations. 

~ C -- i0n - COttpares * evaluates and extends the two approaches for 
helicopter vibration control. 

SecUon_4 shows simulation results based on a simplified model. 

^ L t - i ' ° n -- EiVCS the s ^ a ry. conclusions and proposed work for future 
research. 
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SECTION 2 

nnscRTmoN of tihe-pohatn an d 

FRDlUHNCY-noMAIN fT.TlIOPS 


There are two basic characteristics of the helicopter vibration problem 
that directly impact the active control design and implementation. 


(1) The high-frequency nature of rotorcraft vibration necessitates 
high-bandwidth actuators, sensors and control processors. The 
vibration frequency in most rotorcraft is high enough such that 
aeroelastic and structural modes will no within the control 
bandwidth when active control is used. Aeroelesticity can cause 
rotorcraft dynamics in the neighborhood of the vibration 
frequency to change substantially from one flight condition to 
another. 

(ii) The control action is desired in the neighborhood of widely 
separated discrete frequencies. 

All previous approaches are able to simplify the basic feedback control 
design approach by developing suppression techniques to explicitly control 
disturbances at discrete frequencies. The time-domain approach and the 
frequency- domain approach differ in the manner in which this characteristic 
is utilized, i ,e . , 

(1) The frequency-domain approach attempts to minimize the N/ rev 
(2N/rev, etc.) Fourier transform component of the vibration 
output, through the use of a performance index which depends on 
Fourier components of inputs and outputs. 

(2) The time-domain approach optimizes a performance index with large 
penalty on helicopter response which is narrow-band filtered at 
the vibration frequency. This makes it feasible to use a tine- 
domain dynamic model and update the control law at each sample 
point. 


In this section, the essential characteristics of these approaches are 
described. This description forms the basis of the developments of the next 
section. 
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2.1 FRljtH’I.NO- 1>0 MAIN ROTOIt CRAFT V JURATION CDNTROl.IJiP. 


The frequency-domain approach is baaed on rotorcraft behavior at the 
discrete vibration frequencies. Various formulations can be used. They all 
require a combination of identification and feedback control computation 
steps. 


2.1.1 Model Form 

The control design model for the frequency-domain approach describes 

the relationship between the N/rev harmonics of the control inputs and the 

vibration at desired locations on the rotorcraft. Other dependent variables 

like blade loads may be substituted for vibration. A different model is 

needed for each vibration frequency. A dependent variable nay be a 

performance parameter or any other quantity of interest, c.g., blade loads. 

Using notation of Ref. [3j, let z and 0 be the N/rev harmonics of the 

n c 

dependent-variable vector and the control vector, respectively. In stendv- 
stnte . assuning linearity, z n and O r are related by a transfer matrix, T, 

l „ * T °n + *0 * 

where z^ is the N/rev harmonic of the open-loop vibration level. If the 

linearity does not hold over the entire range of controls, 0 , a locally 

s 

linear aodcl nay be used 


z 

a 


n-1 


T (0 - 

n n 


°n-l> 


( 2 . 2 ) 


In the frequency domain vibration control design approach, vibration 
levels are measured at desired points on the rotorcraft and z is 
computed through harmonic analysis, A three-step procedure is then used (i) 
to identify certain model parameters, (ii) to compute desired feedback gain 
values, and (iii) to compute harmonic components of the desired feedback 
inputs. A schematic diagram of the frequency-domain vibration control lew 
is shown in Figure 2-1. 
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2.1.2 l dent i f lc aJXoji. SJ.e j> 


The identification it a key step in the overall procedure because T 
varies substantially with flight condition, loading end other operating 
variables. z Q depends on the flight condition at well as the gust 
environment and pilot inputs. The model of Equation (2.2) does not 
explicitly take this variation into consideration (it assumes is 

constant). z Q will typically vary core than T. 

The identification problem consists of estimating both T and z^ if 
the model of Equation (2.1) is used. When the model described by Equation 
(2.2) is used, only T needs to be estimated. If T is assumed known or 
does not vary too much, z^ needs to be identified in the global model and 
no parameter needs to be estimated for the local model. 

Various forms of least-squares and Kalman filter may be used to solve 
the identification problem (3). All recursive approaches, including the 
Kalman filter, are of the form (where t is the unknown parameter vector, 

D 

which includes all elements of T and z Q . if nocessary) 


t 

n 


+ K (z - z ) . 
n n n 


(2.3) 


is the best estimate of the parameters based on first n harmonic 

cycles, z is the best estimate of z and can, in general, be written as 
n n 


z -TO ♦ z 
n n n o 


z - z , + T (0 - 0 ,) 

n n-1 n n n-1 


The gain can be obtained by solving estimation error covariance 

equations or be based on least-squares or stochastic gradient equations. A 

Kalman filter formulation has often been used [3]. It is necessary to 

consider T and z^ models driven by white nose such that the gain K^ 

stays finite. Exponential windows may, alternatively, be used to ensure 

that K is not reduced below an acceptable limit, 
n 
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2.1.3 fiain C ompu t a t i on and Control I nput Pc t e mini t i o n 

The io-celled * de tcreii ni st i c controller* computes the feedback gains by 
minimizing a quadratic function of the forn 


z T tf z +0 T ff.O + AO T tr..A0 
nzn n 0 n n AO 


( 2 . 6 ) 


The control law is of the fore 


A0 

n 



where 


(2.7) 


C - -(T T W T + V. + ft' ) _1 <T T W JW.) (2.8) 

z o at* z : v 

assuaing that the first oatriz on the right-hand side is invertible. Note 
T 

that if T W T were invertible, tf. and ff could be set to zero without 

Z AW 

the control activity beconing infinite at any tine point. 

Extensions of the deterninistic controller, described above, include 

stochastic methods, where gains are computed based on errors in parameter 

estimates. This can lead to either a cautious control lew or a dnal 

controller. The main effect of these extensions is to add a positive 

T 

seaidefiuite tern to (T W T + W.„) in Eq. 2.8, thus reducing the gain. In 

z Au 

many cases, a similar effect could be achieved by modifying W or V . . 

t) Au 

(31. 


2.1.4 Implementation of the Erca ue ncy-Doma i n Controller 

The frequency-domain vibration controller involves several steps, which 
are repeated at regular intervals. Figure 2-2 shows the details of the 
implementation for a four-bleded rotor. The input is updated once per 
several vibration cycles (typically once per rev or four vibration cycles 
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for a four-bladed rotor). Tho control update sampling rate is therefore 
once per multiple vibration cycles. 


2.1.5 S ummary 

The frequency-domain approach uses a steady-state model, one of various 
identification procedures (see Section 2.1.2) and one of several gain 
determination procedures (see Section 2.1.3) for the computation of the 
feedback control law. Harmonic analysis methods are needed to determine 
harmonic contents of the vibration at tho N'/rev frequencies. Sine and 
cosine wave reconstruction generates the time domain inputs. Typically a 
control input is applied and the rotor is allowed to reach a pseudo-steady 
state in several vibration cycles. Measurements are then taken for one or 
more cycles of vibration. Vi tlx a small computation delay, the control input 
can then be updated. The control inputs are thus updated every four to ten 
cycles of vibration. Table 2-1 from Ref. [3] summarizes various options in 
the frequency-domain vibration controller implementation. 

Since the control harmonics 0 q are changed continually, the 
rotorcraft may not be in steady state when the measurements are taken. 
Waiting for the rotorcraft to stabilize may cause unacceptable delays in 
computing feedbnek inputs. Thus, an autoregressive moving-average (ARMA) 
form of the model may be more appropriate where z^ depends on past values 
of z and current and past values of 0. Such models are difficult to use 
and have not been studied, though they have been mentioned by Johnson [3]. 
Such models could offer potential advantages in the frequency-domain 
vibration control designs, though they require higher computation time. 


2.2 TIME-DOMAIN APPROACH 

The time-domain approach uses a model defined in the state-space form 
together with a quadratic performance index to design a feedback control 
law. The standard time-domain approach uses a performance index, which 
places equal penalties on states, outputs, and controls at all frequencies. 
A direct application of this approach to the rotorcraft vibration 


9 


suppression control design is not successful because unnecessarily high 
control activity occurs over a wide frequency region. The wide band nature 
of the controller also requires an accurate rotorernft model over a broad 
frer'iency region, a node] which is difficult to obtain because it night span 
many blade and fuselage structural nodes. 

A recent extension of the 1 inear-quadratic-Caussian tine domain 
methodology makes it feasible to place large penalties on outputs at 
selected frequency. This extension is called the frequency-shaping 
methodology [6]. Application of this methodology to rotorcraft vibration 
control leads to acceptable designs. 



2.2.1 Model for the Tine-Oonaii 


The model for the time-domain control design problem starts with 
equations in the state variable form. If u, y and x arc control, outputs 
(or dependent variable) and state time histories respectively, the dynamic 
model takes the form (assuming linearity) 


i •» Fx + Gu , 
y - Dx + Du + w 


(2.9) 


The state variables represent translation, rotational or modal displacements 
and velocities or possibly nonphysical quantities relating inputs and 
outputs, w is the vibration source. The dynamics must include the effects 
of rigid-body and rotor states as well as aeroelesticity and structural 
dynamics states. For a rotorcraft in forward flight, F, G, H, and D 
are periodic functions of time. 

The vibration frequency in rotorcraft is high enough such that the 
blade aeroelastic modes and fuselage structural modes are important. The 
model used for vibration control must include the gain and phase changes 
produced by all the modes up to and beyond the vibration frequency. 

Normally this model could be very complex; however, the narrow-band nature 
of rotorcraft vibration simplifies the model needed for vibration control 
design since the model needs to be accurate near the helicopter vibration 
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frequency. Thu*, a detailed model for each mode is not necessary. A 
narrow-band model could be approximated by a time-invariant reduced-order 
Dodcl. The time-invariant model describes the output in the neighborhood of 
h'/rev resulting from inputs in the neighborhood of N/rev. It would not, 
for example, show harmonics in the output when inputs are applied at N/rev. 





2.2.2 Control Design Approach 

The time-domain approach uses a state-space model and optimizes a cost 
functional that places a large penalty on fuselage accelerations at 
vibration frequencies. The solution to the optimal control problem leads to 
a feedback law where the fuselage accelerations are first filtered by 
undamped, second-order systems. This vibration control solution has been 
possible due to a recent extension of the well-known optimal control 
formulation. The extension allows frequency-dependent penalty functions on 
states and controls [6]. 

Considering the model defined by (2.9), we select a cost functional 
which places large penalties on the vibration output at the vibration 
frequency. The following quadratic cost functional can be used because the 
penalty at the vibration frequency, is infinite. . . 

4 

• w 

(y Ay — 2 V 2 j + u#Bu ) da . (2.10) 

(e> -e> v ) 

where y represents those variables in which the vibration must be reduced. 
Note that there is infinite penalty on the component of y at the vibration 
frequency, y could be components of translational or rotational velocity 
or acceleration at various points on the fuselage. 

The time-domain solution is obtained by defining an additional variable 
vector ( as follows 

2 

yw 

~~i “ 5 ( 2 . 11 ) 

til 

V 
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( 2 . 12 ) 


••2 2 

l 4 w v t ■ “ v y • 

y ' Hx + n , 


and 

J » J U T A£ + u T Bu)dt . (2.13) 

Ihc resulting control law is 

u ** CjX + C 2 5 + , (2.14) 

and the gains are obtained via an optimal control design progran. Because 
of the multiplicity of open-loop eigenvalues, a Schur form (Appendix A and 
Ref. {7}) algorithm must be used. If necessary, the states x may be 
derived from a state estimator based on measurements, y. The method is 
very robust because it car. bo shown that stability will ensure that the 
vibration is completely controlled. 

The basic procedure can be extended to incorporate implementation 
difficulties and other requirements as follows: 

(1) The feedback of rotor states is eliminated by solving for rotor 
states in terms of fuselage accelerations [8]. 

(2) Actuator/ sensor dynamics are included by adding more states. 

(3) The gains can be scheduled with flight conditions, if necessary. 
Extensive scheduling is not likely to be needed because of the 
high margin and xero-to-inf inity gain margin (8). 

Thus, once the model is defined, the entire control design procedure is 
straightforward. The model can be obtained in an online procedure or it can 
be derived off-line. In the off-line procedure, the control gains are 
computed and stored. Only the control input is determined in real time. To 
date, the time-domain approach has only used simulation derived models. 
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2 , 3 Relationship o f the Fr c quencv-Domai n and Time-Domain Mode 1 a 


Tho frequency-domain model of Equations (2.1) or (2.2) and the time- 
domain model of (2.9) are closely related. To demonstrate the relationship, 
we shall show procedures to convert one model form into the other. 

Conversion from the tine-domain to the frequency-donain representation 
is obtained simply by computing the transfer function between y and u at 
the vibration frequency In the following, z n and reprer nt 

complex variables, with real and imaginary parts representing cosine and 
sine components. Note that the time-domain representation results in the 
same transfer matrix for both the cosine and the sine parts. 


x “ y(jw ) ■ 
n J v 


n(jo v i 


F) 


-l 


G u(ju ) + w(jw ) 

. V . v^ 


zo 


(2.15) 


T 0 + Z n 

n n o 


Conversion from the frequency-domain representation to the tine-donnin 
representation is non unique. Without giving any proof, we shall state that 
the simplest tine-domain representation for a general frequency-domain 
transfer matriz is as follows (assuming the order of y is the sane or more 
than the order of u) 


x ■ u , 

y « Hx + Du + w . 


(2.16) 


The order of the state vector equals the number of independent variables 



(2.17) 


Thus, if we choose 
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(2.18) 


it = -w Ib(T ) . 
v n 

D = Re(T ) 
u 

and w(j«D ) = . 

v 0 

The time-domain representation of (2,11) can he iiade equivalent to a general 
frequency-domain model. Note that the open-loop representation is neutrally 
stable. The closed-loop design will, of course, be stable as long as H is 
full rank and all elements of y are to be controlled. 


2.4 SUMMARY 

This chapter summarized the time-domain and the frequency-domain 
methods for active control of rotorcraft vibration. The next chapter 
analyzes each of the techniques to enable an understanding of the relative 
tradeoffs in using the two methods. 
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SECTION 3 

ANALYSIS OK TECHNIQUES 


There are several significant considerations in the use of an active 
feedback control design technique in real systems. These include 
performance under nominal and off-design conditions, ease of design, 
reliability, maintainability, implementation complexity, and robustness. 
Some of the more important analytical considerations are as follows: 


Identification of Mathematical Models for Helicopter Vibration Control 
Desi en - All feedback control design methods require a mathematical 
model to describe the relcti uships between applied control inputs and 
measured outputs. These models may, in general, be developed from 
theoretical analysis or experimental d-ta. For helicopter vibration 
control design, system identifi c ation methods with experimental data 
are desirable because the model is required over a narrow frequency- 
band and the theoretical analysis is likely to be very complex. Both 
real-time and off-line parameter estimation methods must be considered. 

Robustness - The active vibration control system must continue to 
operate satisfactorily and stably with errors in parameters and models 
used for control law design. In addition, the steady-state vibration 
level must be relatively insensitive with respect to weighting 
functions used in control design. 

Transient Behavior - The tine delny between changes in the vibration 
and the rotorcraft approaching steady-state vibration should be small. 
This provides adequate performance in transient flight conditions and - 
in the presence of gusts, which have correlation times of 1 sec. or 
even less. Stability of the control law during the transient must be 
established. As has been shown recently, many adaptive control methods 
may have very poor transient response. 

Susceptibility to Noise - The measurements taken on-board a helicopter 
arc likely to be very noisy because of significant vibration and air 
turbulence. The techniques must be reasonably insensitive to 
measurement and process noise. 


I nnl cmenta 1 1 on Cons i dc r n H o n s - The effect of the active control design 
approach on actuator, sensors, control processor requirements is 
important bocause those will impact the overall cost of the helicopter 
with active vibration control system. 
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During the course of this research, each of these issues was studied. 
In the following sections of this chapter, the most significant results are 
shown. 


3.1 SYSTEM IDENTIFICATION 

This section discusses a selection of algorithms to identify helicopter 
models suitable for use in vibration controller design and implementation. 
Emphasis has been placed on deriving mathematical models from measured 
control inputs and vibration response data. Previous techniques are 
extended to develop low order models that hold primarily in the immediate 
neighborhood of the vibration frequency. Both batch and recursive forms of 
identification methods are discussed. Batch methods are suitable for off- 
line identification, while recursive methods night be applicable to real- 
time identification as well as adaptive implementations. A major portion of 
the discussion is restricted to time-domain model identification since 
identification of frequency-domain models has been covered previously. 

Either an off-line or a real-time (on-line) approach may be used for 
the development of models from test data. If the off-line system 
identification approach is used, three steps are required to derive the 
control law (Figure 3-1). 

In the first step, a test is planned and conducted, where preselected 
multicyclic inputs are applied to the helicopter and the resulting response 
is measured and recorded. Then, the data is used in a batch mode to derive 
a highly accurate model of the helicopter in the operating regions of 
\ interest. The identified model or its simplified form is used to derive the 

control law. Prior to implementation, the control law is evaluated for 
robustness and proper performance over the entire range of operating 

conditions. Failure modes are also tested and redundancy is built to avoid 

i 

catastrophic results in case of failure. Implementation follows, 

‘ In the on-line epproach, all of the above steps are combined. The 

I identification is done on-board in the helicopter and the control law is 

computed using the resulting model. Significant off-line planning and 
analysis are, nevertheless, needed to ensure robustness of algorithms to 
noise,, failures, and sudden changes in model forms. The control law may be 


1 (> 


: j . jc J <*■ - 





updated continuously or at certain intervals. Such methods are referred to 
as * sel f- tuning ' or ’adaptive.* 

Not surprisingly, the basic methods are similar in off-line and on-line 
implementations. In either case, the model needs to be good only in the 
neighborhood of the vibration frequency. The model complexity must be kept 
low both to simplify on-board implementation complexity and because 
unnecessarily complex models can lead to sensitive and nonrobust control 
laws. The key to designing good control laws for helicopter vibration 
reduction is to identify a low-order mathematical model in the neighborhood 
of the vibration frequency and to develop a robust control scheme that can 
use this model. 

The next subsection summarizes system identification methods for the 
derivation of the model for frequency-domain control design. Methods to 
derive applicable time-domain models are discussed in a subsequent 
subsection. 


i 


3.1.1 Identification of Frequency-Domain Models 

Johnson [3] provides a summary of the system identification methods 
that may be used for the frequency-domain models. We shall show the basic 
results for the local model (extensions of the results to other models are 
straightforward) . 


Az = z - z , = TA6 
n n n-1 n 


(3.1) 


Let v be the measurement noise ir z , i.e., 
n n 


z « z + v 
run n n 


(3.2) 


The most common off-line procedure is the least-squares. Substituting 
(3.2) in (3.1), we get 


Az = 7 . - z . 

nm nm n-1 , m 


(z + V ) - (z - V ) 
n n n-ji n-1 
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' / 


= (z - z , ) + ( v - v ,) 
n n-1 n n-1 


of po ;r quality 


«= TA0 + c 
n n 


(3.3) 


where e = v - v .. The estimation model is defined by En. (3.3) with 
n n n-1 

noise e . Note that e is correlated with e , (= v . 

n n n-1 n-1 


V2 ) 


and 


If V is the covariance of the measurement noise v 


z minimizes 
n+1 


C n + 1 (= V n+1 " V- 

and the noise has a Gaussian distribution, the least-squares estimate based 
on • 
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. 

• 

• • 

0 


• 

• 

• 

• 


. 4 V 

-V 


• 

0 

• 

• 


-V 

4V 


e 

l nJ 


(3.4) 


The estimation problem is often simplified by neglecting the off-diagonal 
terms in the covariance matrix. The resulting solution is 



N 


N 

u 

H 2 

Y* (Z - z )V _1 A e T ,. 

n+1, in run n+1 


S 4 Vi v ' 1 49 Ji 


_n=l 


_ n=l 


l-l 


(3.5) 


The most efficient solution is much more complex (see the following). 

The least-squares solution can be converted into a recursive form in 
which the estimate for N+1 measurements can be obtained from the estimate 
with N measurements and estimation error covariance equations. A window 
is usually needed in recursive estimation to ensure that the estimation 
error covariance docs not become too small because a small estimation 
covariance makes the parameter estimates less sensitive to new measurements. 
A Kalman filter can also be used in recursive estimation. For the 
application at hand, the Kalman filter formulation simply formalizes the 
approach to the development of windows in recursive estimation. 
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ORIGINAL PACu 13 
OF POOR QUALITY 

Vc shall Again start with the local model and develop an optimal Kalman 
filter noting the correlation between successive measurements. Ono 
formulation is as follows. 


T 

n 

■ T ♦ 

n-1 

"n 


r 'n 

*= f *= 

V n-1 



* V 

n 



Az 

nm 

c T AO - t, + $ 
n n n 

(3.6) 

n 

n 

is added 

to ensure that the Kalman filter does not 

start 


ignoring future measurements. The formulation of Equation (3.6) is 
interesting for two reasons. First, two additional state vectors are needed 
to model the correlationi secondly, in the extended formulation, there is no 
measurement noise. 

The Kalman filter estimator for the model of Equation (3.6) takes the 

form 


T = T . + (Az - T AO - $ ) K, 
n n-1 run n n n In 


; = ? , + K. (Az - T A9 - L ) 

n n-1 2n na n n n 


t c AZ - T AO -5 
n nm n n n 


(3.7) 


K, is a gain vector (row) and K„ is a scalar. These gains are obtained 
In zn 

by solving the covariance equations. These covariance equations must be 

solved in real-time, in general, because the measurement equations involve 

the input distribution matrix AO . The window size is controlled by 

n 

selecting covariances of noise sources q and v The estimator is 

n n+l 

usually simplified by using the assumption of uncorrelated measurement noise 
in Az 


T = T , + (Az - T A 0 ) K, 
n n-1 nm n n In 


(3.8) 
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I'tt ia«t Io n Error Analysis 

Tho error in the estimate of depends on the Beasurement noise 

covariance assume 1 covariance of » n> and the estimation procedure. 

The error covariance for the exact procedure of Equation (3.7) can be 
determined by solving the covariance equations corresponding to the Kalman 
filter. Usually, the error covariance will be computed for the T matrix 
one row at a tine. 

The estimator model of Equation (3.8) will give larger estimation 
errors than those obtained from the optimal estimation of Eq. (3.7) 

(computed for a single measurement and one row at a tine). In the following, 
it is useful to consider as a row vector: 


T T * T T - K T A0 T T T - r T s ♦ v T, 
n n-1 In n n-1 K ln ? n k ln l n 


(3.9) 


i.V> v. T - h.\ * hn\ 


“ Covariance 
n 


t t 

n 


'n 

L ? nJ 


(3.10) 


Then 


« 0 - 1 'n-1 ’n-1 v n-l 


(3.11) 


where 


(I - K T A0 T ) 
In n 


F T - K T 
In K ln 


(3.12) 
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0 0 E(v v ) 

n n . 


(3.13) 


Equation (3.11) i* solved to determine the error in the estimation of the 

transfer matrix. T, using the simplified Kalman filter of Equation (3.8). 

This estimation error covariance depends on K. A6 and measurement noise. 

Inn 

These error covariances rust be compared with the optimal Kalman filter 
(3.7) to determine if the simplified filter is adequate for vibration 
control . 


3.1.2 Time-Domain Model Identification 


3. 1.2.1 Rotorcraft Model Development 


A general linear model of a helicopter about a trim condition ney be 
written as 


F(v)x + G(,)u 


(3.14) 


where x is the state vector, which may include position, velocity, angles, 
angular rates, rotor states end flexible modes, u is the control input 
vector, and '. is the azimuth location of a reference point on the rotor 
with respect to a reference point in the fixed frame. Also 


(3.15) 


where ft is the rotor speed. For a constant speed rotor, C' » Ot. Hence, 
Equation (3.14) is periodic with respect to time with periodicity Q. 
Outputs y may be written similarly as follows. 
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y e llU’ )x * !)(^)u + w 


(3.16) 


where w it the open-loop uncontrolled vibration. 

A linear helicopter model valid in the neighborhood of the NQ 
vibration frequency can be developed in many different ways (N it the 
number of blades on the rotor). One approach simply computes the Nfl sine 
and cosine components of the output, y, without any input and then with 
sine and cosine inputs at N(l. The data collected in this fora aay be used 
together with the model transformation of Section 2.3 to get an approximate 
time domain representation. N/rev inputs produce higher harmonic responses 
at 2 NO. 3 Nfl . .... etc., which are essentially ignored. This approach 
gives the transfer matrix formulation utilized for much work on higher 
harmonic control [2]. The formulation ignores the transients caused by 
changes in control levels at higher harmonic frequencies. Thus, controls 
based on these models must have update intervals much longer that the 
transient time constant. Johnson [3] proposes a more sophisticated model 
which partially accounts for control transients. 


x 

n 


N 

E 

i"l 


N 

E 


“i Vi + T i Vi 

1=0 


(3.17) 


where z^ represents the sine and cosine components at NQ in the output. 

The time-domain representation is simplified in two stages. First, 
modes which are very far from the vibration frequency are eliminated through 
model reduction. Several methods are available to achieve this model 
reduction. Secondly, if the inputs is of the form u = u^ cos(NQt + <) and 
we are interested in the responses at and around NQ, above equations can 
be simplified to a time invariant model, which well be written as 


x - Fx + Cu + r*(NQt) (3.18) 

y “ Ilx + Du (3.19) 

Note that all matrices in Eq. (3.18) and (3.19) will, in general, be 

functions of In the above derivation, we assume that u and <f> are 

o 

constant or slowly time-varying. 









OF POOrt QUALITY 


The model of Equations (3. IB) and (3.19) can be identified using the 
maximum likelihood approach. The maximum likelihood approach selects model 
parameters F, G, II and D by minimizing the likelihood function, which 
is a function of the innovations sequence and the estimated innovations 
covariance. Hie likelihood function is defined as follows: 



where Vj., k “ 1, 2, .... n is the innovation sequence end B is the 
estimated innovations covariance. The matrix B can be considered as a 
weighting matrix. The innovations are obtained from a state estimating 
Kalman filter (chapter 13, Kef. [9]). 

If the maximum likelihood described above is used for estimation, the 
resulting model will have all the rotorcraft modes, which are excited by the 
inputs and measured by the instruments. Since the model is desired around 
the N/rcv vibration frequency, the maximum likelihood method should be 
extended such that innovations in the region of the vibration frequency are 
given a higher weighting (innovations correspond to differences between the 
estimated model and the measurements). This can be achieved by 
reformulating the likelihood function in the frequency donain by using 


v^(ju) B * V^tju) (3.21) 

max 


Parsevel's theorem. 

u 

max 

L(jw) = 

i ** -u> 


where is the Fourier transform of the innovations and Vj is its 

complex conjugate. To emphasize errors in the neighborhood of the N/ rev 
frequency the likelihood function should be extended to include a frequency 
varying term s(ju), where s(ju) is large at N/rev and small elsewhere. 


L(jtij) * 


u 

max • -i 

£ [V (ju) B V (Ju)]s(Ju) 

“ -U) 

max 


(3.22) 


rs 




m ( jo.) is, thus, a weighting function. Selection of s ( j <u) will determine 
model fidelity (e.g., frequency range over which the model is valid, 
complexity (e.g., model order), and estimation error. An example of s(ju) 

is 


s( jw) 


to 

,2 2 - , .2 
(u - u + .1 UU ) 

V V 


The extended maximum likelihood method can be implemented using the 
following algorithm. 


Step 1: Define a new innovation V as 

Vjtjcu) - Vj( ju) s 1/2 (jw) (3.23) 

Step 2: Develop a Kalman filter whose outputs are V^(jcu). Since the 

model is linear, any frequency shaping applied to u(jw) and 
y(ju) will lead to the sane shaping in the innovations. 

1/2 

Step 3: Pass u and y through a filter represented by s (jw). 

Step 4: Use these nodified n and y for parameter estimation. 

This procedure is shown in Figure 3-2. 


3. 1.2. 2 Model Form and Model Order Selection 

The model form used in identification can be quite general. In 
particular, the specific forms for F, G, T, B and D can be selected 
for any canonical structure, which can be described as a general set of 
linoar input-output relationships. Three of the forms are shown here. Each 
of theso structures represents a state model with minimus number of 
parameters for a given modol order. 
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D is general 



Vutorc gressive Movinp.-Averagc (ARMA) Fora 


If y is the output and u is the input, a discrete ARMA form is 
sinilar to Equation (3.17) except tho outputs are used directly rather than 
the harmonic components at NO. 

Suitable model order depends on the required model fidelity and the 
nature of helicopter dynamics in the neighborhood of N3. Many tests have 
been proposed for selection of model order [10], If the ARMA-type 
representation is used, an adaptive ladder filter might be suitable to 
select model order [11]. Tho baseline approach is to overspecify the model 
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order. The extra poles generally cancel tho extra zeros in the identified 
model . 


Selection of Mod e l From for Time-Do m ain Vibration Control Design 

The autoregressive moving average form and controller canonical form 
have larger universal errors than the real diagonal form. Thus, they are 
suitable when the time-domain rotorcraft model needed for vibration control 
is of low order. The maximum desired order for these forms depends on the 
accuracy of the computer used for system identification and to implement 
vibration control, but will usually be less than ten. The numerical 
conditioning in the real diagonal form does not degrade with increasing 
tiodol order. Thus, this form is preferred when a higher order rotorcraft 
model is desired. 

Simple identification algorithms may be used with controller canonical 
and autoregressive moving-average (ARMA) forms. If significant structural 
modes move around the vibration frequency as the flight condition is varied, 
the model order may have to be estimated in real-time. Simple algorithms 
exist to estimate model order for ARMA forms [11], 

Recent work in ladder-form realizations of ARMA models provides an 
approach to alleviate numerical conditioning problems in ARMA models. 

3. 1.2. 3 Estimation Error Control 

The estimation error will be controlled by selecting a time record over 
which the model for the helicopter does not change appreciably and by 
choosing the filter in Equation (3.23), such that its bandwidth is neither 
too narrow nor too broad. If the filter is too narrow (in the limit a 
single discrete frequency), the level of the signal is decreased and the 
estimation error increases. If the filter is too wide, the signal includes 
dynamics of no interest to the vibration control problem. Even though a 
theoretical analysis can be performed to analyze the effects of filter 
bandwidth on estimation error (see Section 3. 1.2. 4), the selection of the 
bandwidth will vary from one rotorcraft to another. The selection will be 
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guided by the frequencies of structural Bodes, sensor errors, variations in 
rotor speed and rotorcraft transient behavior in gusts and turbulence. 


3. 1.2. 4 Error Analysis 

The error in cstinate of the model order and system parameters is 
determined by postulating a true model, the noise sources and the range of 
frequencies over which the model is desired. 

Let the true model be described in terms of the response of the 
helicopter to sinusoidal inputs at various frequencies. For any flight 
condition, the rotorcraft response amplitude and phase at each frequency 
depends on input amplitude as well os phase because the rotorcraft has a 
nonlinear behavior. Mathematically, the true model is of the form 

y(jw) = T t (je.', u(jw)) + noise (3.26) 

The noise is a combination of sensor errors and gusts. We shall assume that 
a linear model of the following form is desired. 

y (j«o) = T ( jcu,p) u(jai) (3.27) 

ID Z3 

where a parameter vector p will be selected to best fit the measured 
response. We can define a new error term to combine the process and 
measurement noise and the d ifferences in model form 

y( jw) » T (jfc.p) u(ju) + n(jw) (3.28) 

m 

Let N(jio) be the spsetrua of process and measurement noise and modeling 
error. Since we are primarily interested in behavior around the vibration 
frequency, significantly more weighting will be given to the model accuracy 
about that frequency. Lot the weighting matrix be W(jtd). 

The signals y and u are passod through a filter with response 
F(jo>) prior to minimization; therefore, the parameter estimation problem 
involves minimization of the following with respect to p. 



1 

2 


I F( jw)y( ju)-T ( ju,p)F( ju)u(Ju) R F( ju)y(,iu)“T ( ju,p)F( ju)u( ju) du 
J to D 

it) 

(3.29) 


where (.) represents the complex conjugate of (.) and R is a 
frequency independent weighting matrix. For optimization, the derivative of 
the above with respect to p is zero (remember that F(jt>>) is a scalar). 

By substituting for y from Equation (3.28) into Equation (3.29) we get: 

(F ( ju)RF( ju) ) |(T m ( jw.p)-T m ( ju.p)) u(ju) + n(ju)| du = 0 

(3.30) 



Expanding T (ju.p) in a Taylor series 

S3 


3T (jw.f) 

T (ju.p) = T (ju.p) + — Ap (3.31) 

m m _ /v 

0 P 

Vi'e get (dependence on ju is implied) 



A B 


i.e., A A p + /B(ju) n(ju)du = 0 


The covariance of A p is 

E (Ap Ap T ) ° A -1 J B( ju) N(j«) B*( ju)du A -1 (3.33) 

It is possible to show that the estimation error is minimized if F(ju) is 
selected such that (sec Appendix B). 


i.e., R 1/2 F( ju) « lN(ju)] 1/2 


(3.34) 
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The symmetric square root is taken for N(jw). Note that Ft ju>) is a 
scalar, while the NCju) is a matrix. This equality cannot be realized in 
practice unless we consider ono measurement at a time. 


3.2 ROBUSTNESS 

The robustness of both the frequency-domain and the time-domain 
vibration control laws are analyzed. The stability and the performance 
results arc studied first for a fixed gain controller. Stability conditions 
for a controller where parameters are identified in real-time are very 
difficult to analyze and are the subjects of intense research in control 
theory [13]. 

3.2.1 Frequencv - Doncin Controller 

In this analysis, we assume that the correct model for helicopter . 
multicyclic vibration output for time interval n is written as 

z = T0 + z n (3.35) 

n n 0 

where z, 6 and z^ are the vectors of sin and cosine components of 
closed loop vibration, applied input and open-loop vibration. T is the 
transfer matrix. We consider two controllers 


Controller 

1 : 

e = -G z ft , 
n n 0 , n-1 

(3.36) 

Controller 

2: 

A0 » -G z , 
n n n-1 




6 “0 , + A0 , 

n n-1 n-1 

(3.37) 


z_ . represents the estimated open-loop vibration at time n-1 . z , is 

U t fl”l D 1 

the measured closed-loop vioration. The first controller may be referred to 
as the open-loop control law while the second has the closed-loop form. In 
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the following section, stability and noise sensitivity of each of these 
controllers is analyzed. 


Stabi 1 i t v C onditions for Fixed Rain 

In a fixed-gain controller G^ is constant. To analyze the stability 

of the open-loop control law, it is necessary to study variations in Zq 

caused by pilot inputs, gusts, and changes in flight conditions. Let 

z„ , z and z be the true, the estimated and the estimation error in 
Goo 
n n n 

open-loop vibration at cycle n. 




then. 


(3.38) 


,*'o ,»**o -*o- ra .*o -’Vo , (3 - 39> 

n-1 n-1 n n n-1 n-1 


For constant open-loop vibration, z Q 


E(z ) = (I-TG)z. - TG E(z_ ) 

Q U U 

n-1 


(3.40) 


If the bias error in the estimation of z^ is zero, then thv average value 

n 

and variance of the residual vibration are given by the following equations: 


E(z ) = (I-TC)z 
n n 


(3.41) 


Var (z ) = Z « TG EU„ zj )G T T T = TG N , G T T T (3.42) 

n n u . o . n-i 

n-1 n-1 


where N , is the variance of the estimation error in open-loop vibration. 
n-1 

Tho mean value of the vibrntion is small if (I-TG) is small. If T were 
square and invertible, G could be selected to be T * to reduce the mean 
value of vibration to zero. The same would also hold if T had more 
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columns than rows and it were full rank. In addition, the closed-loop 
vibration level is proportional to the estimation error in open-loop 
vibration. 

With a constant G, the open-loop controller will not become unstable 

as long as z is uncorrcl ated with z n _ 1 - A poor selection of G will 

n-1 n 

lead to insufficient vibration reduction or even an increase of vibration. 

Correlation between z^ and z n . can be analyzed as follows. If 

n-1 n ^ 

A represents the projection of z^ on z n _j» then z^ can be 

n-1 n n-1 


written in the form 



Az + noise 
n+1 


(3.43) 


where the additive noise docs not correlate with 
then modifies to 



Equation (3.42) 


Z = TGAZ .AW + additional tern 
n n-1 


(3.44) 


This equation is stable if the eigenvalues of TGA have absolute value less 
than one. Since 7G is of unit order and A must be much smaller than one 
to ensure system stability. 

The basic dynamics of the closed-loop helicopter vibration control law 
is obtained by substituting the control law of Equation (3.37) in the 
dynamic model of Equation (3 35). 


z 

n 


n-1 


TAG = -TG z , 
n n-1 


(3.45) 


If z , is estimated without error, the stability of the closed-loop 
n-1 

dynamic’ depends upon the following 

Abs (eig(I-TG) ) = Abs (l-eig(TG)) £ 1 


The eigenvalues of TG must lie within the unit circle shown in Figure 3-3. 
If the number of outputs to be controlled exceeds the number of available 
inputs, TG is rank deficient end some of its eigenvalues are at the 
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origin. The difference between the number of inputs and outputs is the 
number of zero eigenvalues (Figure 3-3). 

To understand the stability conditions, it is necessary in particular 
to analyze the behavior of the eigenvalues of TG, that are not at tho 
origin. The eigenvalues at origin will stay there because the rank of 
TG can never exceed the rank of t for any value of G. Suppose first 
i that T is known exactly and that the G is selected to minimize the 

vibration to the extent feasible 

G = (T T QT) _1 T T Q (3.46) 

where 0 is a weighting matrix. 

T -IT 
TG ® T(T QT) T Q 

(TG) . (TG) = (T(T T QT) -1 T T Q) (T(T T QT) -1 T T Q) = T(T T QT) _1 T T Q = TG 

(3.47) 

2 

For any matrix where (.) = (.), every eigenvalue is zero or one [12]. 

The right eigenvectors of the eigenvalues at unity lie in the subspace of 
T, e.g., they represent the parts of the vibration which are directly 
controlled from the input. 

(TOT = [T(T T QT) _1 T T G] T = I.T (3.48) 

T 

and the left eigenvectors lie in the subspace of T Q. The right and left 
eigenvectors for unit eigenvalues descrioe errors in T that cause the 
largest perturbation in closed-loop eigenvalues. This will be discussed 
next. Using the small matrix perturbation theory, the perturbation in the 
eigenvalues of T at unity based on perturbation in G is given by [12]. 

(3.491 


(3.50). 
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AX = (T T QT) 1 [T T 0(TAG)T] = AGT 
If T the estimate of T is written as 
T = T-AT 





then 

ACi *= (T T QT - AT T QT - T^OAT + A^QAT)” 1 (T T *AT T )0 - (T T OT)" 1 T T Q 

= (T T QT) -1 AT T a - (T T QT) _1 (AT T UT4T T 0AT)(T 1 QT)T 1 U (3.51) 

Ignoring second and higher order terms in AT . Therefore 

AX *> -(T T QT) _1 T T QAT (3.52) 

The perturbation in eigenvalue therefore depends on error in estimating the 

T T 

transfer matrix as well as the condition of T QT. If T QT is nearly 
singular a small relative change in AT can cause a very large perturbation 
in the eigenvalues. 

One approach to reduce the sensitivity is to add a positive definite 
T 

matrix to (T QT) in the expression for C. This gives the core common 
control gain 

G = (T T QT+A) -1 (T T Q+B) (3.53) 

Now TO has the same number of eigenvalues at the origin since the ranks of 
T and G have not changed. The remaining eigenvalues are not necessarily 
at unity. To ensure that the remaining eigenvalues are at unity A and D are 
related as follows (with this condition TG . TG = TC). 

BT = A (3.54) 

This appears to be a good choice since it places the eigenvalues farthest 
away from the unit circle. The effect of a small perturbation with the 
above condition on the eigenvalues of TG is as follows 

AX ■= (T T QT+A) _1 (T T Q+B)AT (3.55) 

If A is selected properly the condition number of the matrix to be 
inverted can be improved leading to more robustness. 
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Ot her Model ins Er rors 

In the previous analysis, only errors in T were studied and the model 
A7. “ TAO k vas assumed to be the true model. In rotorcraf t there is one 
more error source introduced by pilot and gust inputs. These inputs change 
the uncontrolled vinration level. Let us then assume that the true model is 
of tho form 

Az. » TAO. •* Ax . (3.56) 

k k ok 

Note that Az pk represents the difference between the uncontrolled 
vibration level over two consecutive control updates and can be thought of 
as the rate of change of the uncontrolled vibration. This error can be very 
large if the pilot input or gust spectrum has frequency content in the 
neighborhood of the control update rate, even if the actual variation in 
uncontrolled open-loop vibration is snail. The closed loop dynamics ere 

Az k -TG z k ' :♦ Az ok (3.57) 

z k = (I-TG)z k _ 1 + Az ok (3.58) 

If Az Qk has frequency coaponentt near the poles of (1-TG), then the 
open-loop variation Az Qk can be further amplified in the closed loop 
vibration * k in the steady state. This could seriously degrade the 
performance of the closed-loop control law, if the natural frequencies of 
the external vibration matches the closed loop frequency. 

3.2.2 Time-Domain Controller 

The robustness of the tine domain controller is analyzed moat 
effectively by using robustness theory which is most effectively applied in 
the frequency domain. Let G(jn) and C(ju) represent the transfer 
function equivalents of the rotorcraft and the vibration controller in tbe 
frequency domain (see Figure 3-4). Tbe closed-loop system of Figure 3-4 is 
stable if the encirclement count of the nap det[l ♦ G(ja,) C(j«)J around. 
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the origin, evaluated on the Nsquist contour equal* the number of unstable 
pole * of (iC 1 14 ] . 

Assume now that the controller is stable for the tr-ninal plant. Let 



written as 


C (jtr) - (I ♦ e AG(ju)) G(jui) (3.59) 

where 0 ( t < 1. It is easy to show [13] that the closed-loop system is 
stable if 

o [ G ( j u ) (I ♦ G(jw) C(jb))" 1 ! < 1/ o ( AG( jtu) ) (3.60) 

for all 0 < w < • where o(.) represents the maximum singular value of a 

matrix. The above equation represents a sufficient condition but not one 

that is necessary. Therefore, it may represent a conservative bound. 

Nevertheless, structure of the inequality is a major aid in understanding 

the robustness of the tine domain controller. 

The tine donain vibration controller is designed such that C(ju) has 

a complex pole pair at ♦ jtu^. Therefore, C(jiL') is theoretically infinite 

in the neighborhood of the vibration frequency. It has been observed that 

the optical control formulation selects controller transfer function such 

that G(jiu ) C( jm ) is nearly real and positive. If we assume that G(jte) 
v v 

is a stable transfer function, then because of the frequency- shaping 
selected in the vibration controller design, C(ju) will be small away from 
the vibration frequency. The Dode and Nyquist charts for the vibration 
controller will then be as shown in Figure 3-5. Recause of the phase 
behavior at the vibration frequency and a small closed-loop gain elsewhere, 
the Nyquist chart will be mostly to the right of +1. This shows that the 
phase angle of the rotorcraft transfer function must change by about 90° in 
the neighborhood of the vibration frequency before the controller becomes 
unstable . 
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3.3 TRANSIENT RESPONSE 

The transient response of the vibration control system determines the 
rate at which the initial vibration and the effect of fusts and other 
disturbances arc cliainatcd. A fast transient response is desirable 
consistent with good steady-state performance, stability and robustness. 

This section shows achievable transient response settling times with both 
the frequency-domain and the time-domain controllers. 

The transient settling time will be combination of the time required 
for system identification and for feedback control assuming both are 
performed on-line. The delay in each of the two steps should be isolated to 
understand the weak elements in obtaining fast response. 

Transient requirements can be very stringent. In a gusty environment 
the transient settling tines should be much shorter than the correlation 
times of the gust field. Otherwise vibration caused by gust inputs will 
never be reduced or eliminated. 


3.3.1 F requency- Domain Approach 

The control input in the frequency doaaiD approach is updated after - 
several vibration cycles (typically once per rev or four vibration cycles 
for e four-bladed rotor). The control update sampling rate is therefore 
once per several vibration cycles. Thus, any disturbance which changes 
vibration may take at least one rotor cycle before compensation starts'. 

The transient response of a discrete update system is defined by the 
eigenvalues of the closed-loop dynamics matrix. The slowest component of 
the closed loop response corresponds to the largest eigenvalue of the closed 
loop dynamic matrix (except the eigenvalues on the unit circle). 

Xmax - |Xnax (I-TG)| (3.61) 

If T is the update time, the vibration will be reduced to within 5% of the 
steady state value in 

1 n( .05 ) *t/ ln(Xmax) + average delay in update « -3t/ln(Xnax) + x/2 

3f> 



t/2 represents average delay between the update intervals. 

The key, therefore, is to keep the eigenvalues of (I-TG) as close to 
zero as possible. That requires an accurate knowledge of the transfer 
matrix T because the value of the gain G is computed using the estimated 
value of T. Thus, if T were in error the eigenvalues of (1 -Tj) will be 
nonzero. 


3.3.2 Tine-Domain Approach 


Let the open-loop input-to-output and disturbance-to-output models for 
the helicopter be G(jw) and G^fju). The vibration controller it designed 


with transfer function C(jtc). 

The closed-loop disturbance response is 


G .( I+GC) 
d 


(3.62) 


In the time domain approach the measurements are passed through an nndamped 
filter with frequency at + j« v . If the measurements are y, the output of 
the filter is (see Equation 2.11) 


S + (i) 

v 


and the control law may be written as (see Equation 2.14) 


u - c i* + c 2 ? + c 3 t 


Let the transfer function between the estimated state and the output, y, 
be Tj ( s ) . Then 


(C + C-s) 

u - C lTf (s) ♦ -j [ y 

fa) + fa) 

v J 


CjTjtsHa 2 + w 2 ) + <C 2 ♦ C 3 s) 

2 2 
f + fa) 

V 
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Since f has a pole pair at (1+GC) * hat a zero pair there. The 

closed-loop design places a pole pair approximately at -a ♦ j ui y . Assuming 
that there are no rotorcraft poles in the neighborhood of the vibration 
frequency, the closed-loop behavior for vibration disturbances can be 
approximated to 

k<* 2+ «\> 

G (I+GC) « —= r- (3.63) 

w — >w (s +2$o> s+u ) 

there K is some gain. The response of the closed-loop system to input 

cos (u t) is 
v 

E exp(-£ u t) [cos(u t) + £ sin(o> t)J (3.64) 

V V V 'v V 

Thus, the time constant for the vibration to reduce to within 5% of the 
steady state value is 


U <■> ) 

v v 


(3.65) 


Thus, if u ■ 100 rad sec * and £ «■ 0.1, the vibration will be reduced 

v v 

to within 5% of the steady state value within .3 sec. 

This may seem to imply that the convergence time can be reduced 
arbitrarily by increasing Increasing { y( however, expands the range 

around where the closed-loop gain is large. This will make the system 

less robust with respect to modeling errors far away from the vibration 
frequency. 


3.4 SUSCEPTIBILITY TO NOISE 



The rotorcraft is subject to a variety of disturbances. For the 
purpose of studying the behavior of vibration reducing control laws, all 
external or internal inputs applied to the rotorcraft away from the N/rev 
frequency will be considered as unwanted disturbances. 
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This section describes typical rotorcraft disturbances and then 
analyzes the effects of these disturbances on the behavior of closed-loop 
vibration control systems. 
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3.4.1 P otorcraft Disturbances 

Three classes of rotorcraft disturbance inputs are of interest: 

(a) external gust and turbulence inputs; 

(b) pilot and stability augmentation system (SAS) inputs, and; 

(c) measurement noise sources. 

Ye shall describe each of these disturbance sources to study their effects 
on rotorcraft vibration controller. 

External gust and turbulence arc stochastic inputs in the low frequency 
region of the spectrum. Assuming a Drydcn or VonEarman spectra, the flow 
field variations apply a random force nnd moment input to the rotorcraft. 
Most of the power in the spectra is typically below 1/2 Hz. Gusts nnd 
turbulence affect the flow field around the rotor leading to an increase in 
the open-loop vibration level. The affected flow field may also impact the 
transfer matrix between the N/rev inputs and the associated rotorcraft' 
responses. The major effect is likely to be the change _n the open-loop 
vibration level. 

The pilot inputs are deterministic tico functions whoso spectrum is 
mostly limited to one Hz or less. These inputs can also change the open- 
loop vibration significantly and affect the transfer matrix temporarily. 
Pilot inputs may also produce significant changes in flight condition which 
have more permanent effects on open-loop vibration and transfer matrices. 

The pilot inputs can cause a low-frequency modulation of the cpen-Ioop 
vibration level. Such ucdulation places demanding requirements on the 
vibration control algorithms. 

Mear.urement noise is a combination of random and systematic error 
sources like bias and scale factor errors. 

In the following, wc will analyze the effects of each of these errors 
on closed-loop rotorcraft performance. 
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3.4.2 Frequency-Domain Approach 


Let the average open-loop vibration in the absence of gusts and pilot 

inputs be w.cos{u t+1) where w. is the vibration amplitude and * is an 
0 v 0 

arbitrary phase angle. The effects of pilot, SAS and gust inputs is 
represented as a time-dependent w^(t). The measurement noise is modeled as 
a random additive term, v. The w^tt) has spectrum mostly in the low 
frequency region. 

Consider a discrete implementation, where the vibration frequency is 
biy rad sec 3 and the data is sampled at n points per vibration cycle. 
Thus, the sampling rate is nu v /2n per second and the sampling interval is 
written as 

A *= — (3.66) 

n« 

v 

The ith time domain sample in the measurement is given by 

y i Ew o (Sr) cos ftr + *) +v i (3 - 67) 

The sin and cosine components at the vibration frequency are determined 
as follows (assuming the averaging is done over kth vibration cycle) 

2 n 

* = ~ £ [y. cos (2ni/n) ] 

c. n f — ' i 

k i»=l 




■ » [ c °‘ Is *• (v) 4 £ ■» (“) eo * (“ 4 ♦) 

2 fj ”, (*f)l <3-« 


(3.68) 


All the summations in the above equations are taken with i varying from 
(k-l)n to kn. If there were no error sources, the result would be cos tf. 
The first two terms produce deterministic errors while the last term gives 
random errors. 

To understand the effects of pilot inputs, consider w^ to be a single 
low-frequency sine wave. The effect of a general w^ can be determined by 
superposition 


w (t) = cos(u.t + ij') 
o d 


(3.69) 


where i j> is phase at the initial time. The three terms in equation (3.68) 
are simplified as follows; 


First Torn; 


i=(k-l)n+l 0 “v 11 


i=(k-l)n+l 


” /2niu, \ 

X) cos ( + 'i’J 

l)n+l ' “"v ' 

testes, 

V \ u v \ ' 

/*«, 2n(k-l)o d \ 

* n sin I + — — — l for w. 

\ u o> / d 

\ v v ' 


for 0 ). << <i> 
d v 

(3.70) 


ORIGINAL PA'J.2 V3 
OF POOR QJALlT/. 




Second Term : 


ORIGINAL PfiQZ ‘S 
OF POOR QUAU r; 


kn 

E 

i e (k-l)n+l 


M «• hr ♦ '•) ■ * '") "• hr * ••) 


0 for w. << iu 
d v 


( 3 . 71 ) 


Third Tern : 


This is a random noise terra. Its mean value is assumed zero and the 
covariance is computed as fellows (where r is the covariance of the random 
noise ) 
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E E 


2 2 
n^ cos 
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( 3 . 72 ) 


because the expected value of a cosine wave over its entire cycle is one- 
half. 

The noise in the harmonic cosine component (and similarly in the sine 
component) consists of two parts — a slowly varying function consisting of a 
sum of sine wave and random noise. The random noise increases the root- 
menn-square (RMS) residual vibration in the closed-loop and will impact the 
identification procedures to some extent. The first error term could cause 
major deterioration in identification accuracy because Az^ could occur 
with zero AG^. 

In addition to the additive noise in the measurement of the average 
open-loop vibration, pilot, SAS and disturbance inputs may also cause 
variations in the input to measured response transfer matrix. No analytical 
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models are available to determine the level of the changes in transfer 
matrix caused by a certain size gust input. These effects arc difficult to 
analyze in the frequency domain formulation but could be significant in the 
performance of the overall controller. 



3.4.3 Time-Domain Formulation 

In the time domain formulation the model is 


x *» Fx + Gu 

y ■ Hx + Du + w (3.73) 



where w represents the time history of the open-loop vibration. Sampled 
measurements at n samples per vibration cycle a'.'e 

y = Hx(iA) + Du(iA) + tr(iA) + v, (3.74) 

B i * 

where is measurement noise. The controller is of the form 


x *= <J; x + G y 

°i+l C c i C “i 

(3.75) 

u( iA) = H x + D y 

c c^ c'm^ 

(3.76) 

equations for the system may also be sampled 


x i - x(iA), etc. 

(3.77) 

Vi = * x i + Vi 

(3.78) 

y. *» Hx + Du^ + w 

(3.79) 


In the time-domain formulation, the entire control law and system are 
defined by a set of linear differential or difference equations. Therefore 
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an appropriate Lyapunov equation needs to be solved to compute steady-state 
root-mean-square vibration levels 1 7 J . The steady-state response could be 
determined more accurately by modeling with an appropriate spectrum. 

An approximate analysis can also be performed in the continuous time-domnir. 
representation. Appropriate Lyapunov equations are easy to derive. 

An adaptive controller or a real-time system identification 
implementation could have difficulties here because of the low frequency 
modulation of the open-loop vibration. This problem could be handled by- 
updating the model only after several sample points, or by placing a lower 
confidence on previous open-loop vibration estimate. 


3.5 IMPLEMENTATION CONSIDERATIONS 

Three issues arc iaportant in implementation: (a) actuator 

requirements, (b) sensor requirene - t s, and (c) computational capability. 

Any approach for active vibration reduction requires actuators which can 
produce sufficient deflection of the swashplatc at the N/rev vibratiou 
frequency. Thus, similar actuators are needed for cither approach. The 
only difference is that the commanded input at N/rev is changed in steps in 
the frequency-domain approach and smoothly in the time-domain approach. The 
sensor requirements are also similar in the two approaches. Thus, the major 
difference may arise because of the computation requirements required in the 
two approaches. 

This section discusses specifically the number of computations which 
must be performed in real-time to implement frequency domain and time domain 
control l.ws. To develop a uniform approach to compare the two approaches, 
the following notation will be used 

N = number of blades 

n *» number of samples per vibration cy-cle 
m «=• number of measurement channels 
q = number of independent inputs 


UU 



3.5.1 Freouency-Donain Control l.aw 


Table 3~1 shows the various steps in the frequency-domain control law. 
The number of computations required in an all-digits! implementation are 
also shown. The computations are based on a Kalman filter parameter 
estimator and a continuous control formulation. The number of calculations 
required for other frequency-domain formulations is similar. 

3.5.2 Time-Domain Control Law 


Table 3-2 shows the number of calculations required for a time domain 
controller implemented in discrete form if all identification is performed 
off-line. 

The additional number of calculations for on-line parameter estimation 

is given in Table 3-3 for various approaches which »re applicable for 

helicopter vibration reduction (s is the number of state variables). These 

calculations include the control computation as well. 

The fixed-gain time domain controller requires very few computations 

for implementation. The number of computations required for scheduled gain 

control law is also similar. The computation time will go up by a factor of 
2 

s where s is the number of states. Depending upon the complexity of the 
helicopter model and the possible presence of structural modes in the 
neighborhood of the vibration frequency, the increase could be an order of 
magnitude . 

3.5.3 Comnar i sons 

The time-domain approach requires fewer computations thin the 
frequency-domain approach if a fixed gain (or scheduled gain) controller can 
be used in the time-domain formulation. The real-time computations r.re 
about the same for the time-domain and the frequency-domain controller if 
the model needed for the tine domain control must be identified on-line. 
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SECTION 4 

SIMULATION RESULTS 


This section discusses simulation results. The results have been 
organized to address specific issues in the implementation and use of the 
time-domain and frequency-domain methods. 

Section 4.1 describes the model used in the analysis of this chapter. 
The data is based on a wind-tunnel test performed at NASA 40 z 80 Foot Wind 
Tunnel. Every attempt has been made to obtain results with wide 
applicability. The procedure for comparison and evaluation is summarized 

Section 4.2 covers system identification methods and results achievable 
under noisy conditions. The following section describes the effect of 
closing the control loop assuming model parameters are known. The 
simultaneous identification and control problem is discussed in Section 4.4. 
Section 4.5 gives a short summary. 


4.1 HELICOPTER MODEL AND EVALUATION APPROACH 


The model used in the analysis is based on helicopter vibration data 
collected in the NASA 40 z 80 Foot Wind Tunnel (McCloud and Chopra [15]). 
Table 4-1 shows transfer matrices calculated from ezperimental wind-tunnel 
data for the Raman 7.02-m diameter rotor. It is a four-bladed rotor turning 
at 300 rpm giving the vibration frequency of 20 Hz. Higher harmonic control 
was applied through the controllable twist rotor. The data was collected 
with 2/rev, llrev and 4/rev inputs in the rotating framj. For the purpose of 
our modeling we assume transfer matrices are available at 3/rev, 4/rev and 
5/rev inputs. Note that the transfer matriz coefficients corresponding to 
the 2/rev input and the longitudinal response are very small and unreliable. 

A state variable model was developed using the cosine snd sine 
components cf the response. $, G, H and D matrices for one possible 
choice of the state variable model, which reproduces the transfer matriz, 
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arc shown 1 a Table 4-2. The model was used to simulate responses of the 
helicopter to eiulticyclic inputs. Discrete models arc used throughout. 
Models based on the cosine components were used in design while models based 
on both sine and cosine components were used in validation. 

In the analysis presented here, it is assumed that 4/rev vertical, 
lateral and longitudinal harmonics are measured every 6° of the rotor 
arimuth angle. Thus, 60 data points are collected for every rotor 
revolution or 15 points for every vibration cycle. The effect of 
measurement noise level on estimation accuracy is studied parametrically. 

The frequency-domain and the time-domain approaches are compared on a 
one-to-one basis. The comparison consists of three parts: 


(a) Identification step with no closed-loop control 

(b) Control design step (assuming the model is estimated a priori), 
and 

(c) Closed-loop controller with real-time identification and control. 


The time domain formulation is based on the frequency-domain model 
obtained from the wind tunnel tests. The evaluation procedure, therefore, 
has been extended as shown in Figure 4-1 to accommodate this model form. 


4.2 SYSTEM IDENTIFICATION 

System identification techniques -re used in the time domain and the 
frequency domain. The corresponding results are compared. 

The 4/rev inputs are applied simultaneously in all channels in the 
nonrotating systems. The phase and amplitudes are changed randomly in all 
three channels at the end of each rotor cycle. The maximum amplitude is 
about .005 in the nondimensional units represented by the data (Figure 4-2). 
Data is collected for 25 rotor cycles giving a total of 1501 points over 5 
seconds. Noise free rotorcraft response to this input is shown iu Figure 4- 
3. 

Noise free measurements ere studied first. Then three cases with 
measurement noise and noise due to pilot inputs are evaluated. Thus, four 
cases arc studied in all. 
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(i) hoisc-free measurement* (Figure 4-3) 

(ii) Measurement* with white random noise (noise RMS is 5S. of the 
open-loop vibration level). The response is shown in Figure 
4-4. 

(iii) Measurements with white random noise (HMS is 25* of the open- 
loop vibration level). The corresponding measurements are given 
in Figure 4-5. 

(iv) A 0.85 Hz pilot input or gust disturbance, which causes the 
open-loop uncontrolled vibration to vary * 25* about its nominal 
open-loop value (0.85 Hz is selected to be less than 1 Hz and to 
not coincide with any of the vibration frequency subharmonics). 
The variation in open-loop vibration is shown in Figure 4-6n. 

The response of the rotorcraft is shown in Figure 4-6b. 


4.2.1 Frequen cy-Dona i n Formu l ation 

Cosine and sir.e components are extracted from the measured responses. 
AZj. arc computed fron the sice and cosine responses, h'ote that one out of 
four vibration cycles is used for identification, because the rest of the 
vibration cycles are used to ensure that the output has reached steady state 
prior to measuring the response. 

The correct transfer matrix as well as the estimated values after 25 
cycles arc shown in Table 4-3. 

The rate at which parameters converge for 5* noise are shown in Figure 
4-7 (two sample runs). A corresponding plot for 25* noise level is given in 
Figure 4-8. The convergence tine can be more than one second. 


4.2.2 Time -Dona i n Fornul a ti on 


A batch maximum likelihood technique is used in the time domain. 
Identified parameters for the four cases are given in Table 4-4. Figures 
4-7 and 4-8 shows convergence of parameters in the licc-donain model. 
Convergence time is substantial and will depend upon the shape of the filter 
used. 
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4.3 EVALUATION OF THE CONTROL LAW 


Several control laws ar-. signed both in the time-donai 
frequency domain. The effect measurement noise and pilot 
studied. trade-offs between residual vibration ar.d speed of 
cttabl isbed. 


n end the 
inputs is 
response are 


4.3.1 Frcquencv-Poce i n Formulation 

The closed-loop systen is sinulated by collecting data over one 
vibration cycle, setting aside one vibration cycle for computation and then 
applying the computed feedback input. Nothing is done for the following two 
vibration cycles. Thus tbe feedback control signal is conputed and updated 
once every rotor cycle. 

In the following, two sets of plots are shown for esch case. The first 
plot shows the sensor measurements, which include both the measurement noise 
and the rotorcraft vibration. The second plot shows the actual vibration of 
the rotorcraft. While the second plot shows the vibration tbe crew and the 
airframe experiences, the first plot indicates the quality of signal used 
for control. 

Sine and cosine components z. of the vibration along the three axes 

i 

are entracted. Based on the previous measured values of Az^ * s 

computed. This is multiplied by the gain matrices to compute inputs for 
feedback . 

The following gain matrices are used. 
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(1) G * 1 , since there are three inputs and three controls (each 
with a sine end a cosine component) in our eiample, it is possible 
to make this choice. This is the fastest possible controller, 
leading to vibration reduction in about four vibration cycles, 
under ideal conditions. 

(2) G * nT where a < 1. Tho time constant increases as a is 
decreased from one. 

(3) G *= (t'oT ^ A) * (T T Q * li) , 0 is set to identity, and A and B 

are varied. 

The results for measurements with and without noise are shown in 
Figures 4-9 to 4-19 and arc summarized in Table 4-5. As a general 
conclusions as the gains are reduced the steady-state vibration response due 
to measurement noise goes down hut the convergence time goes up. The 
presence of pilot-induced open-loop viDration can significantly degrade the 
closed-loop results and high gains are needed to minimize the impact of this 
effect on closed-loop performance. 




4-3.2 Timc-Docair. Fo rm ulat ion 

In the time domain formulation, the basic dynamic model of the 
rotorcraft is extended to include frequency-shaping of helicopter vibration. 
This gives a 9-state formulation for the three-axis problem. The control- 
law is designed by optimizing a quadratic cost functional in the states of 
the intended system and the input. The model in the discrete formulation 
has the form: 
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The nominal value of h i% 3 z 30 . The control law it of the form 


u *■ C x ■* C x ' ♦ C x * 
i 1 i Vi 3 i-1 


Table 4-6 shows the closed-loop eigenvalues of the tine domain 
controller with b * 3 » 10 , The absolute value of the closed-loop 
eigenvalues is also shown. As is well known for digital systems, the 
distance of the closed-loop eigenvalues away from unit circle is a measure 
of system convergence time (the farther the eigenvalue the better the 
convergence time). The convergence tine for the open-loop vibration to 
reduce to 5% for the three coaplex eigenvalue pairs is as follows: 


-3/ ln( .989) . -3/In(.9766) and -3/ln(.9737) sanple points 

or 771. 123 and 113 saaple points. Thus, the worst case convergence 

time is 271 sample points. 

Vith sampling tine of 300 per second, the convergence tine for various nodes 
varies between .37 5 and .98. To this must be added the state estimation 
delay to determine the overall convergence tine. 

Figure 4-19 shows the response of the closed-loop system when there is 
no measurement noise. The convergence time varies from 0.3s to 0.6s. 

Figures 4-20 and 4-21 show the helicopter vibration response and measured 
outputs in the presence of 5“» white Gaussian noise. Note that the 
measurement is very noisy even in steady state (representing mostly 
measurement noise). The vibration is reduced to a lew value within about 
the same convergence time as before. 

Corresponding rotorcraft vibration and measured time histories for 253> 
measurement noise are shown in Figures 4~22 and 4-23. Note that there is a 
significant increase in the closed-loop vibration level due to the noise. 

The response of the closed-loop system is computed with a pilot, SAS or 
gust input which causes a modulation of the open-loop rotorcraft vibration. 
The modulation has a magnitude equal to 25^> of the open-loop vibration with 
a frequency of 0.85 IIz (see Figure 4-6s). Figure 4-24 shows the closed-loop 
response of the system. The response is modulated at 1.70 IIz and has a 
magnitude of about 10T of the open-loop value. 
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The doted- loop responses for the baseline controller arc summarized in 
Tabic 4-7. The RMS value of residual vibration for noisy measurements is 
also determined statistically. actual simulated values are close to the 
statistically determined quantities. For this controller 25% white Gaussian 
measurement noise docs not appear to cause major problems. Pil ot- induced 
responses could be important. 

To study the behavior of the closed-loop system with variations in 

transient response time, two additional control laws are designed. In the 

7 

first one, the control penalty is reduced by a factor of ten to 3 * 10 and 
in the second one by a factor of one hundred to 3 x 10^. These penalty 
reductions vill cause the closed-loop eigenvalues to migrate away from the 
unit circle causing the transient response to be faster. Figure 4-26 shows 

7 

the closed-loop response with 5% noise level and b * 3 r 10 (medium speed 
controller). The response of the medium speed controller with pil ot- induce d 
variations in open-loop vibration is shown in Figure 4-26. A comparison of 
Figures 4-20 and 4-25 shows that increasing the speed (or gain) of the 
controller increases the residual closed-loop vibration due to measurement 
noise. A similar comparison of Figure 4-23 and 4-26 indicates that the 
reverse is true of pilot induced variations ir. open-loop vibration. Thus, 
in the design of closed-loop control system for rotorcraft vibration control 
a trade-off must be established to ensure that reasonable residual vibration 
is obtained both in the presence of measurement noise and open-loop 
variations in rotorcraft vibrations. 

The magnitude of the eigenvalues and the corresponding convergence 
times in term.' of the number of sample points is given in Table 4-8. 
Transient response can be significantly accelerated by simply decreasing the 
control weighting. 

Table 4-9 shows the RMS steady-state vibration level for the three 
controllers with 25% measurement noise and with pilot-induced variation in 
open-loop vibration. The faster response causes an amplification of the 
measurement noise. It is very beneficial for pilot induced input. 

The controllers with faster response are also less robust. Modeling 
errors farther away from the vibration frequency can drive the controller 
unstable or degrade performance. 
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4.4 REAL-TIME ESTIMATION AND CONTROL 

This section discusses the impj (mentation of vibration control laws 
where the systeia identification and feedback control functions are performed 
in real-time rather than off-line. The issues of significance are 

(i) Stability problems in the presence of noise, 

(ii) The intervals over which the model must be updated. 

(iii) Steady-state response, and 

(iv) Need for caution or dual crntrol. 

These areas are studied through simulation in the following subsections. 


4.4.1 Frequency-Domain Formulation 

The frequency-domain control law is simulated using the Salman filter 
estimates and gains from Section 4.3.1. Results in the presence of three 
different kinds of noise sources (discussed in Section 4.2) are shown in 
Figures 4-27 to 4-29. Note significant performance degradation when pilot- 
induced inputs modulate open-loop helicopter vibration levels. 


4.4.2 Time-Domain Formulation 

The time-domain approach is implemented as follows: 

• Control law updated every rotor revolution and is based on the 
solution to a Riccati equation 

O Identification performed recursively, updating parameter values at 
every sample point 

•Estimation accuracy computed using information matrix 
•The entire procedure is implemented in square-root form 

• Sampl ed-da ta control design and covariance analysis 
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Results in tho absence of measurement noise are shown in Figure 4-30. 

Figure 4-31 expands the first 500 data points and shows the building of 
control and filtered state. 

A wliite-Caussian noise with an RMS value of .6 is added to the 
measurements. The measured time histories are shown in Figure 4-32. The 
actual accelerations with closed-loop control are given in Figure 4-33. 
Closed-loop system appears well controlled. 

The corresponding time histories for noise within + 1/rev around the 
vibration frequency (a two pole filter is used with white noise input) are 
shown in Figure 4-34 and 4-35. Note that there is significant helicopter 
vibration in the steady-state. Part of the steady-state vibration is due to 
parameter estimation error and partly because of feeding measurement noise 
directly into the plant. 
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SUCTION 5 

SUMMARY. CONCLUSIONS A N D FUTURE WORK' 


The analysis presented in this report compares the frequency-domain and 
the time-domain approaches for effectiveness, robustness and implementation 
complexity. The study has attempted to point out the advantages and 
disadvantages of both approaches. 


5.1 CONCLUSIONS 

The following conclusions may be drawn based on the theoretical 
analysis presented here and the simulation results based on one specific 
model form. 


System Identification Methods 

System identification methods are available for both the frequency- 
domain and the time-domain approaches. These techniques may be implemented 
on-line in a recursive mode or off-line in a batch mode. The off-line 
approach can give more accurate models for specific flight conditions. To 
track parameter variations, the on-line approach will have to use a past- 
fading approach. 

It is necessary to base models on dots in the immediate neighborhood of 
the vibration frequency. Appropriate approaches which make this possible 
have been specified for the time-domain approach. 

Vhen local models are used in the frequency-domain approach the pseudo 
measurement noise is non-white. Use of appropriate noise models can give 
improved results. 

t-rrrr:o p.\r.r RI.AXK not filmed . ■ 
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Robustness 


The time-domain approach appears to be more robust (has higher gain and 
phase margins). Most of this robustness is available mainly because the 
inputs can be updated more frequently. However, it appears that gain 
scheduling can be used in the frequency-domain approach as ir. the time- 
domain approach. More points may be needed in the frequency-domain approach 
than in the time-domain. 


Transient Response 

When the rotorcraft is operating at one-flight condition and the open- 
loop vibration is constant, transient response in the frequency-domain is 
superior to that in the time-domain if the transfer matrix is known exactly. 
The settling time in the frequency-domain is about one rotor cycle. If any 
of the above conditions do not hold, the transient settling tine can 
increase substantially. The faster settling time also places certain 
closed-loop eigenvalues nearer the unit circle with possible impact on 
robustness. The transient response for the tine-domain controller may be 
made arbitrarily fast. However, faster response reduces robustness because 
it increases gain awny from the vibration frequency. In practice, the 
transient settling time will be determined by a priori confidence in the 
model and the distance of structural modes from the vibration frequency. 

Susceptibility to Noise 

Roth techniques are highly susceptible to noise, particularly 
variations in open-loop vibration due to pilot inputs and noise. The 
residual vibration increases as the transient settling time is reduced. 
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I mp! enenta t i on Comrlciit v 

Hoth approaches require similar actuators and sensors. The time-domain 
controller is easier to implement even with on-line identification. 


5.2 FUTURE RESEARCH 

Future research should be performed along the following directions: 

(i) Establish definitive relationships between robustness, transient 
response and susceptibility to noise for each of the two 
techniques. 

(ii) Study the effects of other kinds of noise sources. 

(iii) Further simulations for both approaches. 

(iv) Perform wind-tunnel and flight tests for each of tin two 
approaches to provide an understanding of the impact of 
unmodeled parts of the system on closed-loop behavior. 
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TABLE 2-1. SUMMARY OF PAST WORK ON REGULATORS FOR HELICOPTER VIBRATION ALLEVIATION 
(Taken from Ref. 13]. All references may be found in that report.) 


Peculator 

type 

Reference 

Investigation 

Identification 

Controller 


Invariable 

Kretz et al. (1973a, b) 

Expel imont 

beast squares 

J = z T z 


open loop 

McCloud and Kretz (1974) 

Experiment 

Least squares 




McCloud (1975) 

Theory 

Least squares 

J = z T V,' 7 z 



McCloud and Weisbrich (1978) 

Experiment 

Least squares 

J » z T W ? z 



Brown and McCloud (1980) 

Experiment 

Least squares 

.1 = «■ ■■V..v 



Sissingh and Dunham (1974) 

Experiment 

Direct Inverse 

Direct inverse 



Powers (1973) 

Experiment 

Various methods 

Direct inverse 



Wood et al. (1980) 

Experiment 

Various methods 

Direct inverse 


Invariable 

Shaw and Albion (1980) 

Experiment 

Direct inverse 

Direct inverse 

/ * a * *•* y 

closed loop 

Shaw (1980) 

Theory 

Direct inverse 

Direct inverse 


Adaptive 

Hammond ( I 980 ) 

Experiment 

Kalman filter, T and z 5 

.! » ••• T W.z, caution 

* t r \ r , a - 

Open loop 

Molusis ct al. (1981) 

Experiment. 

Kalman filter, T and z e ; 

.1 " z^W ? r, caution 

Rate I;-:*.. 




or only z 0 



Adaptive 

Shaw (1930) 

Theory 

Kalman filter, local T 

Direct inverse 


closed loop 

Tavlor, farrar. and Miao (1980) 


Kalman filler, local T 


:t « l-C rev. 


Taylor ct al. (I960) 

Theory 



rate I i“it . * 


"Included in the theoretical development, but not iscd in applications. 
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TABLE 3-1. COMPUTATIONS REQUIRED TO IMPLEMENT FREQUENCY DOMAIN CONTROLLER 

(LOCAL MODEL APPROACH) 


Step No. 

Description 

Relevant Equations 

No. of Computations* 
Per Rotor Cycle 

1 

Harmonic Analysis 

V 2at \ 

z ic iL/ cos n 1 Equivalent 

t=1 I to FFT for 

and ( n fe w 

nN , . ( spectral 

z. » V y, sin — — 1 l incs 

is *— ’ •'i n I 

L=1 1 

i “ 1 , 2, 3 . . . m 

2nmN 

(Note: FFT will 

require more compu- 
tations, but is not 
necessary) 

2 

Feedback Controller 

°k ' Vk-l + C 2°k-l 

4qm 

3 

Control Time 
History Generation 

cos 2ut/n 

u t " 

sin 27tt/n 

2nqN 

U 

Identification 
(Simplified Kalman 
Filter Formulation) 

? k ’ Vi + w vW K I 
"k ' "k°k /( V 40 X°k> 

"k ' »-Vi 10 I.i>Vi«-Vi m I.i> T 
+ "i-AiViVi 


5 

Cain Computation 
(Cautious Con- 
trol ler) 

. 

C » -dt t w 
1 z 

C 2 « dwao 

D - <t\t + W Q + W A0 )' 1 



One computation is approximately one multiplication phis one addition. 
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TAIlI.i: 3-3. NIIMllKR 01- COMHIJ ATIONS lt»R Hil.I.Y-AliAI'I'l VI- 
TIM-.- DOMAIN (Y)NTKOl. LAW 


\ 



Id update 

Computational 
Compl ex i ty 

Numerical lichavior 

Suitability for adaptive 
control implementation 

Unmodified Kalman 
update 

<6s 2 *8s+l)m 

poor, results in 
'burst instability* 

Numerically unacceptable 

II— 1> factored 
covariance update 
with regularization 

(6s 2 +8s)ft 

information 

good, inappropriate 
for fixed-point 
ar i thac tic 

Easy implementation of 
minimum variance tyre i 

control laws, difficult to ; 

1 

exploit a-priori 




Can use a-priori structure 
and information directly. 
Can be used with variety of 
control design algorithms. 

I’ottcr aquare root 

(12s 2 +7 s+3)n 

reasonable 

needs square roots 

Good 
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TAHI.I: 4-2. STAll SI’ACi: Dl-SCKll'riON MATKK1S LOR 
ROTOIUHAIT V1RRAT10N AT 80 KNOTS 


! - 0.4000 0.0000 o.oooo" 

0.0000 0.4000 0.0000 

0.0000 0.0000 0.4000 

IUted on 

111 

104.1444 -134.8268 40.4742 

-97.5373 -52.4049 26.9036 

1.0457 -98.4880-198.2593 


112 

" 95.3033 -147.2091 39.4285* 

-38.9294 -79.3323 -52.5713 

2.9470 -305.4933 -195.6688 


(I - 0.6000 0.0000 0.0000 

“ 0.0000 0.6000 0.0000 

0.0000 0.0000 0.6000 


ni 

173.7178 -81.1690 52.7490 

-56.9288 -45.0650 4.3788 

4.6005 -25.6997 -196.6914 


1)2 

132.1210 -90.7718 53.0735 

-24.5194 -80.6299 -43.9430 

3.7123 -165.9383 -151.6213 


cos 

stn 


0.1634 

0.6845 

Vertical 

-0.4094 

0.2145 

Lateral 

0.0269 . 

0.0362 

Longitudinal 


where 


s. . * . + 

G .u. 


l+l i 

d i 


>i “ H, x k 

+ D„ u. : 

1 i 

cosine components 

>i = »2 «i 

+ n 2 u. : 

sine components 
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TAMP 4 3. 1 PINT 111 ID VAMUS 01 TKANS1TK MATRICI-S 1-tlR VARIOUS CONDITIONS 


4-3a. No Measurement Noise (Simulation Values) 


248.4905 
177.9708 
81 .8083 
-59.2212 
76.6687 
-23.0155 


-126.9578 

-82.6902 

23.6948 

55.4642 

29.7998 

-15.2986 


5.3513 

-96.4113 

-339.0359 

-0.5946 

56.0047 

112.7393 


59.2212 
-76.6687 
23.0155 
248.4905 
-177.9708 
81 .8083 


-55.4642 

-29.7998 

15.2986 

-126.9578 

-82.6902 

23.6948 


0.5946 

-56.0047 

-112.7393 

5.3513 

-96.4113 

-339.0359 


4~3b. 5% Measurement Noise 


228.7618 

-201.2055 

86.2400 

-34.8128 

49.0249 

-38.7759 


-115.5355 

-79.3397 

19.3265 

64.6834 

18.5928 

-23.0295 


4.1620 

-96.4760 

-337.7200 

- 1.0120 

56.6752 

113.7915 


61.1591 

-62.2130 

20.5434 

239.5251 

-176.1633 

67.2966 


-79.8408 

-16.3239 

18.8274 

-123.9286 

-84.5502 

35.6074 


-0.0975 

-55.7717 

-114.1909 

4.2645 

-97.9689 

-339.7685 


4-3c. 253b Measurement Noise 


318.8766 

-226.8184 

48.7592 

-139.5581 

21.1616 

-27.7391 


-102.2032 

-51.7225 

29.2625 

20.8387 

82.6690 

-36.9317 


3.6299 

-101.0769 

-334.6284 

-1.0341 

57.3932 

117.2422 


94.2875 

-91.5136 

76.3503 

153.3878 

-150.7690 

73.8549 


-77.8795 

-12.6879 

-12.6542 

-138.1573 

-60.6836 

51.1235 


3.0921 

-56.5297 

-109.4967 

8.5748 

-96.8024 

-339.7697 


4-3d. Pilot-Induced Variation in Oprn-Loop Vibration and No Noise 


243.2912 

-159.9086 

87.1926 

-69.9638 

62.5147 

-26.4943 


-146. 3046 
-10.3621 
44.2526 
18.9536 
-18.9634 
-24.5629 


5.7786 

-98.0988 

-339.4992 

0.1511 

57.0134 

112.8793 


88.4342 

-186.9508 

-8.1353 

302.8968 

-105.1557 

95.0282 


-54.3660 

-33.0340 

14.2206 

-124.2955 

-79.2573 

24.8475 


2.4589 
-63.0788 
-114.7309 
8.7989 
-91 .7920 
-338.2183 



















TAiu.r 4-5. n.osi:i»- i-ooi* itki-ormanci- or various ikloulncy* domain 
VIBRATION CONTROL LAWS UNUI-R VARIOUS C0NHJT10NS 
( VA1.ULS AH!- OIVLN FOR LAC1I OF 1111- HlKl-K ONl-S) 


Vibration 

No Noise 

. 

5% Measurement 
Nui sc 

25% Measurement 
Noi se 

1*11 ot- Induced 
Variation in 
Open-Loop 
Vibrations 


Sensor 
Out rut 

Vibration 

1 

Sensor ' 

Output Vibration 

» 

Vibration 


- 

.0373 

.0123 

.1867 

.0615 

.0631 

G = T _1 

- 

.0246 

.0072 

.1232 

.0358 

.0424 

l 

J 


.0023 

.0009 

.0116 

.0046 

.0042 

i 

.0132 

.0379 

.0145 

.1796 

.0394 

.0840 

G *= 0.5T _1 

.0090 

.0247 

.0086 

.1183 

.0184 

.0559 


.0008 

.0024 

.0010 

.0112 

.0030 

.0053 

G = (T T T+A) _1 







(T T +B) 

Unstable 






O 

rH 

a 

< 







j B = 400 

t . _ 







f 

j G = (T T T+A) _1 

.0107 

.0406 

.0202 




<T T +B) 

.0322 

.0704 

.0668 

Not 

Not 

Not 





Con- 

Com- 

Com- 





pated 

puted 

ruted 

A = 105 







B = 0 

.0012 

.0047 

.0042 






TA1II.H 4-6 KlGI-NVAI.Ul-S 01- nil- CLOSF.IV-I.OOI* SYSTFM 
Willi b - 3*10 1 ' 


Absolute Value 
of Closed- Loop 


Ope n- 

■Loop 

Closed- 

•Loop 

Eigenvalue 

0.9135 

4 

0.40671 

0.9035 

♦ 

0.4025 i 

0.9890 

0.9135 

4 

0.40671 

0.9035 

- 

0 .4025 i 

0.9890 

0.9135 

4 

0.4067i 

0.8918 

4 

0 .3980 i 

0.9766 

0.4000 

4 

O.OOOOi 

0.8918 

- 

0.3980i 

0.9766 

0.4000 

- 

O.OOOOi 

0.8891 

+ 

0.3970i 

0.9737 

0.4000 

+ 

O.OOOOi 

0.8891 


0.3970i 

0.9737 

0.9135 

- 

0.4067 i 

0.4000 

4 

O.OOOOi 

0.4000 

0.9135 


0.4067 i 

0.3999 

4 

O.OOOOi 

0.3999 

0.9135 

- 

0.4067 i 

0.3999 

4 

O.OOOOi 

0.3999 
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TABLE 4-7. CLOSED- LOOP VIBRATION 



Open- 

Loop 

Vibration 

5% Measurement Noise 

23% Mens 

urenent Noise 

With Pilot (or Gust) 
Induced Variation 
in Open-1. oop Vihratio* 


Based on 
Sample 

Statistical Iv 
Cal ciliated 

Based on 
Samp 1 c 

Statistical ly 
Calculated 

Sample 

Based 

First Output 

.7037 

.0^49 

.0053 

mom 

.0265 

.0575 

Second Output 

.4622 

.0015 

.0025 

.0074 

.0125 

.0496 

Third Output 

.0431 

.0004 

.0004 

.0019 

.0018 

.0042 

















TAHl.li 4 8. COMPARISON OI- Till- CONTKOLLl-RS V111I 71IRICK 1)1 IHCRKNT 
HCNALTIICS ON CONTROL INHIT 


. Kigenval ue 
Magnitude 


Convergence 
Time to 5*fc 
of Open-Loop 
Vibration 
(sample points) 


High Control 
Weighting 
b * 3 x 10 

Medium Control 
Weighting 
b = 3 x 10 7 

.9890 

.9659 

.9766 

.9292 

.9737 

.9210 

1 271 

86 

123 

41 

113 

36 


Low Control 
Weighting, 
b - 3 x 10 


f 







TAlil.E 4-9. STIiAHY- STATE RMS VI li RATION FOR CONTROU.F.RS IN 1111-. 
PRESENCE OF MEASUREMENT MUSE AND PILOT- INDUCED VARIATIONS IN 

OPEN-LOOP V lit RATION 


Disturbance 

Outputs 

Slow Controller 
b = 3 x 10 

Medium Speed 
b -- 3 i 10' 

25$> White 

First 

.0265 


Gaussian 

Second 

.0125 


Measurement 

Third 

.0018 


Noise 



RH 

Pilot Induced 

First 

.0575 


Variation in 

Second 

.0496 


Open-Loop 

Third 

.0042 

.0019 

Vibration 
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Figure 2-1. A Schematic Diagram of a Frequency-Domain Vibration Controller 
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Figure 4-Jlb. Actual Rotorcreft Vibration Using Frequency-Dona in 
Control with 25% Measurencnt Noise and G c T * 
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Figure 4-12. Rotorcraft Vibration with Frequency-Domain Control Law, 
Pilot-Induced Variation in Open-Loop Vibration and G = T _1 
(Note: Each second corresponds to 5 rotor cycles) 
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Figure 4-14*. Frequency-Domain Controller with 5% Measurement Noise 
and G “ 0.5 T * (Sensor Outputs) 

(Note: Each second corresponds to 5 rotor cycles) 
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Figure 4 _ 14b. Continued (Actual Vibrations) 
(Note: Each second corresponds to 5 rotor cycles) 
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Figure 4-lSa. Frequency-Domain Vibration Control Law with 25% 
Measurement Noise and C = 0.5 T _1 (Sensor Outputs) 
(Note: Each second corresponds to 5 rotor cycles) 
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Figure 4-15b. Continued (Actual Vibrations) 
(Note: Each second corresponds to 5 rotor cycles) 
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Figure 4-16. Actual Rotorcraft Vibration with Frequency-Domain 
Controller in the Presence of Pilot-Induced Variation in 
Open-Loop Vibration and G = 0.5T 1 
(Note: Each second corresponds to 5 rotor cycles) 
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Figure 4-17. Frequency-Domain Vibration Control Law with 5% Measurement 
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(Note: Each second corresponds to 5 rotor cycles) 
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Figure 4-20. Closed-Loop Response with Measurement Noise 
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Figure 4-24, Closed-Loop Response of the System with Pilot Inputs 
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Figure 4-27 . Measured Rotorcraft Vibration with Joint Estimation and 
Frequency-Domain Control (5% Measurement Noise) 
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Figure 4-28. Measured Rotorcrcft Vibration with Joint Estimation and 
Frequency- Domain Control (25% Measurement Noise) 
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Figaro 4-29. Measured Rotorcraft \ibration with Joint Estimation and 
Frequency-Domain Control (Variation in Open-Loop Vibration Due to 

Filot Applied Inputs) 
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Figure 4-35. Closed-Loop Vibration with 
Narrow-Band Heasuresent Noise 
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APPENDIX A 

MATRIX": A DATA ANALYSIS, ^ySTFJt _IPF.NTI£IJC&ILQ^ 

CONTROL DESIGN. AND SIMULATION PACKAGE 

MATRIX^ is an interactive software systen for the conputer-aided-design 
and analysis of control systems for dynamic plants. In particular, data 
analysis, system identification, model redaction, control design, evaluation 
and simulation steps can be performed conveniently using MATRIX^. A user- 
friendly interpreter incorporates powerful operations and easy-to-use 
graphics. The systen can reduce control systen design costs and provide 
capability for rapid analytical, laboratory and full scale testing of new 
control concepts. 

DESIGN PHILOSOPHY 

The design and development of MATRIX^ is based on two sinple ideas: 

1. It should be possible to nso nost design and analysis techniques to 
solve any problem. The systen must be responsive to the particular 
tastes of a control systen designer. The outputs should be 
available in formats that are easy to conprehond. 

2. A single computer program nust be capable of perforning nost 
operations of interest to a control design engineer, because 
transferring files and data from one program to another can consume 
a najor portion of an analyst's tine. 

MATRIX^ implements commonly used design and analysis methods and 
facilitates using other current techniques and even algorithms which night 
be developed in the future. The software doos not advocate any particular 
philosophy for control design or systen identification. For ezanple, both 
transfor-function and state variable control system design techniques are 
offered and can be applied to tho sane problem with oqual ease. The overall 
flexibility of the package will become apparent as the software is 
described. 
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(•RJiCiilMXt; PACE plant; not filmed 


MATRIX^ provide* a comprehensive tet of capabilities in a tingle 
integrated package with unifora data structures and file fornats (tee Figure 
A-l). Most 'bookkeeping* chores are perforated by the software, leaving the 
control designer free to tackle control problcns. It is even possible to 
implenent linear control lews and estinators designed by MATRIX^ directly 
into a multivariable control processor, MCP-100, designed and built by 
Integrated Systens, Inc. This iaplenentation requires no knowledge of 
assembly language or other real-tine progressing techniques. 


MATRIX STRUCTURE 

The structure of the software has been optimized for flexibility, ease 
of use, expansion capability and to reduce the need for the control systems 
experts to learn the details of operating systems and programming languages. 
Command inputs and graphics or alpha-numeric outputs are in clearly 
understandable fornats. 

Figure A-2 shows the overall structure of MATRIX^. The program uses a 
stack architecture for storing variables in current use, allowing excellent 
interactive response. Command files and data files, not in current use, are 
stored on disc. Macros and current data are available on the stack. Filo 
management is completely user transparent. 

Macro and command file capabilities are used to build a hierarchical 
algorithm structure in MATRIX^. The cost basic (and often the most used) 
algorithms are programmed as language primitives. Primitives are the 
easiest to use and will oxecute rapidly. Examples of primitives are 


Macro capability allows a particular string to be replaced by another 
(usually much shorter) string. This capability must exist in all good 
interactive programs or operating systems. 
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arithmetic operations on matrices, singular value decompositions, Ricceti 
equation solution, fast Fourier transform (FFT) and square root estinate 
update. Specialized and loss conmonly used operations are written as macros 
or command files. When research leads to new algorithms, these algorithms 
are conveniently written as macros and command files during the 
'experimental ' stago. Only a selected set of algorithms should then be 
added as language primitives, if desired. A command file or a macro 
consists of a series of language primitives and other Qacros and command 
files, 

MATRIX^ borrows certain aspects of its architecture from MATLAB, a 
matrix laboratory developed by Moler at Stanford University and now at the 
University of New Mexico, Albuquerque. MATLAB is in the public domain. 


MATRIX CAPABILITIES 

A 

The overall capabilities of MATRIX^ may be divided into the following 
broad categories: 

1. Matrix, vector and scalar operations, 

2. Graphics, 

3. Control design, 

4. System identification and signal processing, and 

5. Simulation and evaluation. 

The details on how each of these capabilities is utilized are shown in 
subsequent sections. A brief summary is given here. 

Matrix, Vector and Scalar Operations (Table A-l ) 

Basic arithmetic operations on compatible matrices are performed using 
standard symbols for addition, subtraction, multiplication, division end 
raising to a power. For example, the square root of a matrix A is 
computed by tbo command 
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SQRT(A). 

Basic algorithms for linear system solutions, eigensystca decomposition 
(including reliable determination of the Jordan form), singular value 
decomposition (SVD), QZ decomposition, and matrix algebraic operations are 
implemented as language primitives. Many of the primitives are inherited 
from MATLAB. Most commonly used operations on matrices are available. 


Graphics (Table A-2 ) 

A powerful, flexible, user-driven graphics capability is available in 
MATRIX^. The command, PLOT(x,y), causes y to be plotted against x 
with automatic selection of scales, axis labels and area of the screen where 
the plot is to be made. Only a part of the screen is used for plotting, 
leaving the remainder for alpha-numeric output and command inputs. 

Using English language commands, it is possible to change the size of 
the pints, x-axis or y-axis scales, location of the plot on the screen, axis 
labels, title, grid lines, tic marks and other variables. The available 
command structure makes graphics capability extremely useful and user- 
friendly. 

Control Design (Table A~3) 

In MATRIX^, control design can be based on any of the following: 

(a) Classical methods including Bode and Nyquist (including methods 
for multivariable plants), 

(b) Linear-Quadratic-Gaussian (LOG) approach, 

(c) Methods based on A-3 invariant subspaces, 

(d) Eigenstructuro assignment and zero placement, end 

(e) Adaptive control using self-tuning regulators and other 
techniques. 
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For the LQG problem, the algebraic Riccati equation la aolved from 
extended Hamilton equations avoiding inverses, which are troublesome in the 
singular case. The equations are row compressed with an orthogonal 
transformation followed by the QZ pencil decomposition and a backward stable 
ordering of the eigenvalues. 

Meaningful extensions to LQG methods require inclusion of dynamics of 
reference, disturbances, sensors and actuators. Appending of dynamics in 
frequency-shaped control design or model-following techniques involves 
forming augmented equations, which are easily accomplished with MATRIX^ 
primitives. Use of frequency-shaped costs with singular value plots for 
robustness evaluation allow incorporation of engineering judgment in the 
control design. 

Evaluation tools for linear systems include frequency response, power 
spectral density plots, time responses, transmission zeros and individual 
transfer function zeros. The principle vector algorithm (PVA) primitive for 
numerically reliable extraction of the Jordan Form (with discriminatory rank 
deflation of root clusters) is very useful in modal analysis of open-loop 
systems of vehicles and structures. PVA permits computation of residues or 
partial fraction expansions of multivariable systems. 


Data Analysis and System Identification (Table A~4) 

Data analysis and identification can be performed very efficiently and . 
easily in MATRIX^. Tied with a flexible graphics package, MATRIX^ provides 
a production environment for batch and recursive identification methods. A 
universal interface to external simulations is provided to facilitate data 
transfer. Data can be cencored, detrended and analyzed in MATRIX^. Batch 
procedures include the standard regression methods with analysis of variance 
and step-wise regression. State-space and nonlinear batch maximum 
likelihood procedures are also available. Recursive algorithms such as the 
recursive least squares, recursive maximum likelihood and extended Kalman 
filter with Ljung's modification are available. All covariance 
factorizations and updates aro in D-D form for numerical reliability. Non-, 
parametric batch and semi-batch methods using the FFT are provided for 
auto/cross covariances/spoctras. Adaptive control algorithms for 
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multivariable systems using U-D updates can be designed using simple 
commands . 

Filter design facilities include Finite Impulse Response design in 
Chebyshcv Norn, Wiener Filters, window-based designs, and Infinite impulse 
Response design. 


Simulation and Evaluation (Table A-S) 

MATRIX^ provides capabilities for efficient linear and nonlinear 
simulation. The linear simulation is performed with a discrete 
representation and is structured to fully utilize any sparseness in system 
matrices . 

A general nonlinear simulation model may be directly connected to 
MATRIX^. The software is capable of solving explicit differential equations 
as veil as implicit equations of the form g(x,x,t) “ 0. 

Use of 'chopped arithmetic,' i.e., using various effective machine 
word-lengths, can provide performance evaluation of on-board small word- 
length control system implementations. Simulations with finite word-length 
processing elements can thus be performed efficiently. 


NUMERICAL AI-GORITUHS 

Numerical reliability and stability are important in all of the 
analysis environments described above. Primitive matrix operations are 
based on the best available numerical software drawn from EISPACK, UNPACK 
and recent research in numerical techniques . 

, Reporting and control of numerical errors are comprehensive. All 

f potential loss of accuracy that may be significant is reported to the user. 

Backward stable unitary transformations are performed on data whenever 
possible. 



AVAILABILITY 

MATRIX it currently operational on VAX with VMS and UNIX operating 
systems, IBM 3033 and 3081 with MVS/TSO operating system and CDC series 
machines with NOS operating system. Implementation on other machines is 
in progress. 


now 
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Figure A—l. MATRIX^: A Data Analysis, System Identification, 

Control Design and Evaluation Package 
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TABLE A-l. KATOIX x CAPABILITIES: MATRIX ARITHMETIC 


Data Entering, Display end Editing 
Addition, Subtraction, Multiplication and Division 
Absolute Value, Real Part, Imaginary Part and the 
Conplcx Conjugate of a Matrix 
Sum end Product of Matrix Elements 
El enent-by-El ement Multiply and Divide 
SIN, COS, ATAN, SORT, LOG. EXP of Matrix Elenents 
Eigenvalue, Singular Value, Principal Values, Schur, 
LU, Cholcsky and QR Decomposition of Matrices 
Random Vector end Matrix Generation end Manipulation 


TABLE A-2. MATE I X^ CAPABILITIES: GRAPHICS 


Flexible Commands 

Multiple Plots 

Axis Labels and Plot Title 

Symbols and Lines 

Tics and Grids 

Log Scales 

Bar Charts 

Plot Location and Size 
Personalizing the Plot Command 
Report Quality Plots 
3-D Graphics 

Parallel and Perspective Projections 

Surfaces 

Curves 

Viewing Transformations 





TABLE A- 3 . MATRIX CAPABILITIES: CONTROL DESIGN AND 

SYSTEM ANAI.YS1S CAPABILITIES 


£ii*i4c^_Xpois (applicable to nui tivariable eysteas) 

Bode Plots 
Nyquist Plots 

Root Locus (with Zoon Capability) 

Modern Tool s 

Optical Control Design. Discrete and Continuous 
Optinal Filter Design, Discrete and Continuous 
Frequency-Shaping LQG Design 

Singular-Value Dcconposition of the Return-Difference 
Eigensystes Decortpositions Including the Jordan Cannonical 
ForB 

Model Following Control 
Model Reduction 

Linearization of Nonlinear Systeas 
Minical Realization and Hainan Deconposition 
Geoaetric Control Algorithas 
Multivariable Nyquist Plots 



TABLE A-4. MATRIX^ CAPABILITIES: SYSTEM IDENTIFICATION, SIGNAL 

PROCESSING AND DATA ANALYSIS CAPABILITIES 


Tine-History Plot* 

Multichannel Crocs-Plott 
Scatter Plot* 

Frequency Plot* 

Ilistograas 

Data Transforaations and Spectral Analysis 

Detrending 

Censoring 

Digital Filtering 

Discrete Fourier Transfora 

Inverse Fourier Transfora 

Autocorrelation 

Cross Correlation 

Antospectrua 

Cross Spectrua 

Decication and Interpolation 

Maxines Entropy Spcctrua Estination 

System Identification 

Step-Wise Regression and Model Dnilding 
Maxinun Likelihood Identification of State-Space Models 
and Nonlinear Models (generated by Systen-Build) 

Recursive Maxiraun Likelihood Identification 
Extended Kalaan Filter Algoritha 

Fil ter Design 

Window-Based Methods 
Wiener Filter 

REMEZ Exchange Algoritha for Finite Ispnlse Response Filters 
Elliptic, Chebyshcv, Buttervorth Infinite Iapnlse Response 
Design 



/ 

/ 


TAtiLE A-5. MATRIX CAPAniLITIES : SIMULATION AND ^SYSTEM 

fiuiLDING CAPABILITIES 


Graphical Display of Block Diagrams and Interactive 
System-Building 

Continaons and Discrete Systems 
Linear and Nonlinear Coaponcnts 
Vector Signals and Multivariable Connections 
Variety of Input Options for Linoar Systems 
State-Space 
Transfer Functions 
Pole-Zeros 
Matrix Fractions 

Multirate and Mixed System Analysis 


Interactive system-building will be released in June 1983. 
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APPENDIX B 

OPTIMAL SELECTION OK FREQUENCY SHAPING FOR ESTIMATION 


We want to minimize the determinant of tbe estimation error by a choice 
of the shaping function, P. 


E < A p A p> * A*" 1 J' B( ju)N( jw)B*( jw)du A' 

where 


(B.l ) 



(B.2) 


(B.3) 


Computation of the deaired fora for F(ju) when N(jw) ia a general natrix 
leada to complex equations which cannot be used in practice. We drive the 
results for the case where N( jw) is scaled by the seme frequency function 
g(ju), i.e., all measureuent noises have the same spectrum. Then R will 
be selected such that 


H(jo>) - g( ju)R 


(B.4 ) 




(Fdg( ju)F) du 


(B.5) 


Since A and B are both quadratic functions of F, a constant 
change in F at all frequencies will not affect the estimation error. 
Thus, we can impose the constraint 




(F*F)dw 


(B.6) 


i’RECKM.Vc; I'ACI’ I'.l.W'K N 


°T I’li. M Fry 
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0HIG1NAL 
OF POOR QUALITY 

Since the error covariance la a matrix, we nuat convert it into a 
scalar fora for ainiaizatlon. It ia typically useful to oininize tho 
de terninant , Using Lagrange aultipliers for the above constraint of 
Equation (B.6), we get 

L * det (X)) + XI / F Fdw-1] 

“ 1 2 Idet .f("air) F ' RF ("air) + x J*F*Fde>-l) 

(de t (A) ) 2 ' a P ' J 

(B.7) 

£ 

Taking derivatives with respect to F F, we get 



+ X - 0 


The above equation is satisfied if 

F g(j<«))F “1 at all frequencies 


(B.O) 


(B.9) 




or 


F » [g(ju)] 


- 1/2 


This equality is physically neaningful because it inplies a filter which 
attentuates the aeasureaents over that range of frequencies where 
neasurenent noise is large. 

For a general N(jcj), the above equation is empirically eztended to 
tho following for F g(ju)F - 1, A - C. Thus, the quantity inside tho 
bracket is zero for all w. 


R 1/2 F(j<i>) w [M(js»)j“ 1/2 
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